Magnetic properties of nanoscale compass-Heisenberg planar clusters 
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We study a model of spins 1/2 on a square lattice, generalizing the quantum compass model via 
the addition of perturbing Heisenberg interactions between nearest neighbors, and investigate its 
phase diagram and magnetic excitations. This model has motivations both from the field of strongly 
correlated systems with orbital degeneracy and from that of solid-state based devices proposed for 
quantum computing. We find that the high degeneracy of ground states of the compass model is 
fragile and changes into twofold degenerate ground states for any finite amplitude of Heisenberg 
coupling. By computing the spin structure factors of finite clusters with Lanczos diagonalization, 
we evidence a rich variety of phases characterized by Z2 symmetry, that are either ferromagnetic, 
columnar, or of Neel type, and analyze the effects of quantum fluctuations on phase boundaries. 
In the ordered phases the anisotropy of compass interactions leads to a finite excitation gap to 
spin waves. We show that for small nanoscale clusters with large anisotropy gap the lowest excita- 
tions are column-flip excitations that emerge due to Heisenberg perturbations from the degenerate 
ground states of the compass model. We derive an effective one-dimensional XYZ model which 
faithfully reproduces the exact structure of these excited states. The low energy column-flip or 
compass-type excitations are robust against decoherence processes and are therefore well designed 
for storing information in quantum computing. We also point out that the dipolar interactions 
between nitrogen-vacancy centers forming a rectangular lattice in a diamond matrix may permit a 
solid-state realization of the anisotropic compass-Heisenberg model. 

PACS numbers: 75.10.Jm, 03.65.Ud, 05.30.Rt, 64.70.Tg 



I. INTRODUCTION 



Frustrated quantum magnetism belongs to the very ac- 
tive research areas in condensed matter theory. Frustra- 
tion is one of the simplest concepts in physics with far 
reaching consequences.—^ It is well known that antifer- 
romagnctic (AF) exchange for three spins on a triangle 
is geometrically frustrated, both in classical Ising and 
in quantum Heisenberg models. On a two-dimensional 
(2D) square lattice, frustration typically involves inter- 
actions between further neighbors competing with those 
between nearest neighbors, but it can also occur with 
only the latter ones: for instance when, compared to the 
Ising model, the sign of spin exchange along every sec- 
ond column is reversed. The resulting model, called fully 
frustrated Ising models is exactly solvable, with a phase 
transition at a lower temperature than that of the 2D 
Ising model^ — the low temperature phase, with exten- 
sive entropy due to frustration, is described in terms of 
dimer coverings of the dual latticed A quantum analog of 
this classical frustrated model on a square lattice is the 
2D quantum compass model (QCM)^ where interactions 
couple either Sf or S? components of nearest neighbor 
5 = 1/2 spins, depending on the spatial bond direction 
x or z respectively. When the associated exchange con- 
stants [see Fig.QJb)] have different values, these two spin 
components are nonequivalent. Yet even otherwise, and 
in spite of its (2 + 1) dimensionality, this model displays 
a finite-temperature phase transition of the 2D Ising uni- 
versality class £ but the symmetry broken phase at low 
temperature is characterized by high ground state degen- 



eracy in the thermodynamic limits 

One can interpolate between the 2D QCM and Ising 
models by modifying continuously the spin components 
coupled on the different bonds <2 This allows to highlight 
that the QCM is closely related to orbital physics, where 
the exchange interactions are directional^ In fact, for 
e g orbitals one finds a 2D superexchange model where 
the interactions are gradually modified from the classi- 
cal Ising model toward the QCM. While the frustration 
of interactions is clearly weaker in the e g orbital model 
than in the QCM, the latter may be considered as a 
generic description of frustrated directional orbital in- 
teractions which arise for strongly correlated electrons in 
transition metal oxides with partly filled degenerate 3d 
orbitals. In these systems the orbital degrees of free- 
dom play a crucial role in determining ground states 
with coexisting magnetic and orbital order, described 
by spin-orbital superexchangeJi~— The orbital interac- 
tions are intrinsically frustrated^ as they have low sym- 
metry in pseudospin space representing the orbital de- 
grees of freedom. Typically the symmetry is that of 
the lattice due to the shape of 3d wave functions, and 
not the SU(2) symmetry typical for exchange spin in- 
teractions. Although frustration is at its maximum in 
three-dimensional (3D) models and it was concluded from 
the high-temperature expansion that a phase transition 
to the symmetry-broken states does not occur^ recent 
Monte-Carlo simulations have shown that symmetry- 
broken phases are stabilized at low temperatures both 
in e g and t^g orbital models i 18 ! 19 This contrast with a 
formal analog of the QCM defined on the honeycomb 
lattice, the Kitaev model^ for which an exact solution 
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evidenced a spin liquid ground state. In fact, both mod- 
els can describe orbitally degenerate Mott insulators in a 
limit of strong spin-orbit coupling 2 ^ which selects at each 
transition metal ion a low-energy doublet represented by 
a pseudospin-1/2 variable — which model is actually rel- 
evant depends on the geometry of the system. However, 
realistic orbital models are more involved,— inter alia 
due to non-conservation of the orbital quantum numbers 
that follows both from hybridization processes with oxy- 
gen orbitals in an oxide and from the structure of charge 
excitations controlled by Hund's exchange in orbital de- 
generate systems. The QCM was designed to avoid all 
these complications and to address a paradigm of intrin- 
sic frustration due to directional conflicting interactions. 

Another motivation for introducing and investigat- 
ing the 2D QCM comes from the field of quantum 
computation.— Several realizations of computing devices 
with protected qubits have been proposed, with descrip- 
tions related to the QCM, and in various contexts: (i) 
in Josephson junction arrays , 23 > 24 as well as (ii) with po- 
lar molecules, or (iii) with systems of trapped ions in 
optical lattices £5. In general, in order to construct a de- 
vice which could serve for information storage, a mani- 
fold of degenerate states is required,— and these degen- 
eracies should be stable against noise and other small 
perturbations^ thanks to particular symmetries of the 
Hamiltonian. This is actually the case in the quantum 
compass model, where two types of symmetry opera- 
tions described by operators P\ and Qj commute with 
the Hamiltonian but not with each other, see Sec. II. As 
a result the eigenstates of the system are characterized 
by related integrals of motion, and concerning the ground 
state, by a hidden dimer-dimer order j21 more importantly 
for quantum information, an exact twofold degeneracy 
of all quantum levels was evidenced on finite clusters of 
arbitrary size.— These degeneracies are, thanks to the 
non-local nature of operators P, and Qj , robust to local 
perturbations; in consequence qubits defined by a real- 
ization of the QCM are expected to be protected against 
noise, so that this model is of prime interest for quantum 
computing. 

The QCM has thus an interdisciplinary character as it 
plays an important role in the modeling not only of corre- 
lated transition metal oxides, but also of protected qubits 
for quantum computations. An intriguing question im- 
portant in all these contexts and asked shortly after the 
QCM was introduced concerns the nature of a quantum 
phase transition (QPT) that occurs when anisotropic in- 
teractions are varied through the isotropic point, also 
called compass point. A first order transition between 
two distinct phases with directional ordering, along ei- 
ther rows or columns, was suggested by Lanczos diago- 
nalization and Green's Function Monte-Carlo simulations 
for finite clusters^ and later confirmed using a projected 
entangled-pair state algorithm.— At this transition, a dis- 
crete symmetry in spin space is spontaneously broken 
since the frustrated interactions along two different di- 
rections are equivalent and the spin orientation follows 



one of them. In terms of broken symmetries this transi- 
tion is remarkably similar to the first order QPT found in 
the exact solution of the one-dimensional (ID) QCM,— 
or a compass ladder^ where two different types of or- 
der stem from the invariant subspaces of the ID model, 
which suggests that a similar mechanism may operate 
also in two dimensions. 

Particularly in the context of proposed realizations of 
quantum computing devices based on finite spins clus- 
ters with compass-like interactions, a fundamental ques- 
tion to ask is how the highly degenerate ground states^ 
are modified when a small perturbation occurs. We ar- 
gue that Heisenberg interactions between nearest neigh- 
bor spins stand for a class of perturbations to the com- 
pass terms which are typical in solid state systems - 
for instance a Hamiltonian with compass and Heisenberg 
terms would describe exchange processes in some Mott 
insulators with strong spin-orbit coupling and 180-degrec 
bonds.— In a broader perspective we study in this work 
the effect of Heisenberg perturbations, by considering 
a generalization of the QCM called compass- Heisenberg 
(CH) model. We find that the high degeneracy of ground 
states in the thermodynamic limit (TL) is removed by 
Heisenberg terms of arbitrarily small amplitude, and var- 
ious magnetically ordered phases arise, with a preferred 
spin direction related to the ordered pattern. In macro- 
scopic systems the lowest excitations are thus gapped 
spin waves; on nanoclusters however, for small enough 
Heisenberg amplitude another type of excitations can be 
of lower energy than spin waves: these are the column-flip 
excitations, from the ordered ground states selected by 
small Heisenberg terms to the many other eigenstates of 
the low-energy manifold minimizing the energy of domi- 
nating compass interactions. The column-flip excitations 
arc robust with respect to decay into spin waves, and 
have an original multiplet structure which can be cap- 
tured by an adapted effective model; this analysis leads 
us to propose that in a regime of moderate Heisenberg 
interactions these excitations could be used in a novel 
type of solid-state-based quantum computing scheme. 

Recent experimental developments on arrays of 
nitrogen-vacancy (NV) centers, constituting point-like 
defects in a diamond matrix ) 31 i 32 may bring a further 
motivation to the study of a model with compass and 
Heisenberg interactions. Indeed these defects can be ef- 
fectively described by quantum spins Sij& coupled (un- 
der some conditions) predominantly by the dipolar inter- 
actions of the form: 

H *v = E ^{Si-Sj-SS^S]), (l.l) 

where the 7 axis in spin space is along the spatial di- 
rection connecting spins i and j. These interactions 
are long-ranged, but rapidly decaying with distance; if 
defects sit on sites of a rectangular cluster, the (domi- 
nant) interactions between nearest neighbors are a sum 
of Hciscnberg-like and compass-like terms' Beyond the 
nature of couplings, an aspect which must be taken into 
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FIG. 1: (Color online) (a) Structure of a NV center in a di- 
amond matrix, which can be described by an effective spin; 
(b) four NV centers forming a rectangle are controlled by 
dipolar interactions between the associated spins Si. Dipole 
interaction H^ip involves compass-type terms of various am- 
plitudes J x , J z , depending on the pair considered, and Heisen- 
berg interactions, (c) Ground state of anisotropic compass- 
Heisenberg model represented by the spins of NV centers and 
arranged in a rectangular array, (d) Different low-energy ex- 
citations of this system: a column flip (CF) reverses the spins 
of a whole column, while a spin-flip (SF) reverses a single 
spin, contributing to spin-wave excitations. 



consideration in this context is a possible splitting of en- 
ergy levels for a single NV center; but one can a priori 
consider a situation where these splitting are small be- 
fore the typical energy scale of dipolar couplings; alterna- 
tively, other possible realizations of dipolar-coupled spin 
arrays are conceivable (with e.g. I n = 1/2 nuclear spins, 
in a layered crystal with an orthorhombic unit cell and 
negligible effects of hyperfinc coupling to electron spins). 
We will show that such systems can have lowest excita- 
tions consisting of reversing columns of spins, and which 
could be used for encoding protected qubits, see Fig. [1] 

The paper is organized as follows. In Sec. [II] we intro- 
duce the CH model and state the problem of frustrated 
interactions and possible QPTs. There are two variants: 
the ferromagnetic (FM) and the antiferromagnetic (AF) 
CH model. We first focus on the AF CH model, and 
present selected data for the spin structure factors in 
Sec. IIIII and show that long-range order is induced by 
arbitrarily small Heiscnbcrg interactions. The full phase 
diagram of the CH model is presented in Sec. lIVI — there, 
we provide evidence that some phase transitions occur 
for the same interaction parameters as in the classical 
CH model, while others transition lines arc affected by 
quantum fluctuations, see Sees. II VB H and IIVB 2l Next 
we analyze the features of the FM CH model and discuss 
its phase diagram in Sec. IIV CI Spin wave excitations are 
derived and discussed for different phases in Sec. [V] We 
turn then to the analysis of the lowest energy states of 



finite clusters and show in Sec. I VII that: (i) the ground 
state and the low energy excitations are very well de- 
scribed by an effective ID model which captures the es- 
sential parameter dependence of columnar (or column- 
flip) excitations characteristic of the compass-Heisenberg 
model; and (ii) there exists a parameter range where the 
column-flip excitations are the lowest energy excitations 
and cannot decay into spin waves. The paper is summa- 
rized in Sec. I VIII where open issues and possible exten- 
sions of this work, aiming at experimental relevance, are 
also discussed. 



II. COMPASS-HEISENBERG MODEL 

We consider a model of spins S = 1/2 on the square 
lattice, with axes in the ac plane, labeled here x and z 
after the interacting spin components in the QCM. The 
nearest neighbor interactions are of two types: (i) frus- 
trated compass interactions of amplitudes J x and J z , and 
(ii) Heisenberg interactions with an exchange /. While 
the Heisenberg interaction is isotropic in spin space and 
bond-independent the compass interactions depend on 
the bond direction. On bonds along the a;-axis the x- 
componcnts of spins are coupled by terms af j/jf - +1 (we 
label the sites in a 2D cluster by two indices {i,j}), 
and on bonds along the z-axis the coupling concerns the 
^-components, being of the form o-fjCT* +1 j. For conve- 
nience we use here Pauli matrices a? = {o-p,af,ap} with 
r = such that cr5 = ±1 in the a z basis. 

The CH Hamiltonian reads 'M- 



Here the sums over i and j run over the intervals [l,£ z ] 
and [1, L x ] consistent with either periodic boundary con- 
ditions (PBC) or open boundary conditions (OBC). We 
consider rectangular clusters with N = L x x L z spins 
and L Xl L z < 6 for both PBC and OBC. Another type 
of clusters (considered only with PBC) are clusters tilted 
by it /A w.r.t. the previous ones and containing N = 2L 2 
spins (for L = 3,4), they also accommodate all symme- 
tries of the model. 

The structure of eigenstates of T-L depends only on 
the relative amplitude of parameters J z , J x and / in 
Eq. (|2.ip . Thus, the total space of interaction param- 
eters may be characterized by a point on the spherical 
surface, see Fig. O The compass interactions are then 
described by the related interaction strength, 



J c = y/J* + J§, (2.2) 
and the angle <f> determined by the exchange constants, 
J z = J c cos4>, J x = J c sm4> , (2.3) 



FIG. 2: (Color online) Parametrization of the compass- 
Heisenberg model Eq. (|2.1[) by two angles {</!>, A point in 
parameter space is indicated by the (blue) dot on the sphere; 
interaction parameters /, J x , and J z correspond to its carte- 
sian coordinates on the respective axes, while angles 8 and <j> 
are its spherical coordinates. 

and are represented by a point in the (J X ,J Z ) plane, sec 
Fig. [3] Using this parametrization the Heisenberg inter- 
action I may be expressed by the angle 9: 

tan 8=—. (2.4) 

Jc 

In the following we shall denote by antiferromagnetic 
compass-Heisenberg (AF CH) model the case where J z > 
0, and by ferromagnetic compass-Heisenberg (FM CH) 
model the case J z < 0. 

We introduce some non-local operators playing a cen- 
tral role in the QCM (and, as we will see, in the CH 
model), and defined cither on rows i or on columns j of 
these clusters, by: 

/' IK, ( 2 - 5 ) 

j 

Qj H< (2-6) 

i 

Here Qj rotates all spins in the column j from up- to 
down-orientation or vice versa — this operation is called 
column flip (CF), see FigQJd). Similarly, Pj rotates a 
whole row of spins pointing along ±x direction into ^fx 
direction in spin-space. In the QCM (7 = 0) these opera- 
tors are known&2£ to commute with the compass Hamil- 
tonian Hi-o (i.e., [Pi,%i=o] = 0, [Qj , %i=o] = 0) and 

— restricting now to the first type of (untilted) clusters 

— they anticommute {{Pi,Qj} = 0), accounting for an 
exact twofold degeneracy. 

In presence of Heisenberg interactions, the commuta- 
tors above become non-zero. To evaluate those, it is use- 
ful to consider separately two terms oc I, complementary 
to each other in Eq. (|2.ip : the first term Wj acts on hor- 
izontal bonds, while the second term 'H'j acts on vertical 
bonds. It is now straightforward to show that Qj com- 
mutes with the columnar interaction, i.e., [Qj,Hj] = 



FIG. 3: (Color online) Parametrization of the compass model 
by either couplings J z and J x or by the radius J c and angle 
4> defined by Eqs. (|2.2|) and (|2.3|) . respectively. On this cir- 
cle, isotropic compass points \J X \ = \Jz\ are shown by filled 
(blue) square and diamond for the model with AF J z > in- 
teractions and by empty (purple) square and diamond for the 
model with FM J 2 < 0, while Ising points are indicated by 
dots (J x = 0) or vertical bars (J 2 = 0). The points <j> — 37r/20 
and (j) = 23-7r/20, frequently considered hereafter, are also in- 
dicated by radial (red) bars. 

but not with the Heisenberg interactions on the rows 
(h-ll; 2/0, ]T (2.7) 

i,ae{y,z} 

Hence the column j is here coupled to the left (j : — 1) and 
the right (j + 1) column by a z a z and a v a y Heisenberg 
interactions. 

These operators are related to a formalism which al- 
lows to understand the high ground state degeneracy of 
the QCM in the TL, and which we will describe briefly. 
Here, we consider a situation with strong anisotropy 
\J Z \ 3> \J x \i which calls for a perturbative treatment of J x 
couplings. The unperturbed Hamiltonian contains only 
the J z compass couplings, which select, on a (L Xl L z ) 
rectangular cluster, 2 Lx columnar states. These states, 
where in each column all spins point along the same 
axis z — cither ferromagnetically or antifcrromagncti- 
cally aligned depending on the sign of J z - can be la- 
beled using pseudospin variables f), with 1 < j ; < L x . 
For a given columnar state, the column j will be de- 
scribed by an eigenstate of the pseudospin operator rj, 
cither | tj) = +r/| ti) or I ij) = -Tj\ ij), depend- 
ing on whether the spin at a reference site has the 
orientation +z or — z respectively; orientations of other 
spins in the column follow from its FM (J z < 0) or AF 
( J z > 0) long-range ordered nature. In both cases, the 
operators r? and rj flip all spins of the column j with 
amplitudes given by the respective Pauli matrices. The 
operator t? has actually the same action on a column j as 
the operator Qj = J\ i afj , with the only difference that 
Tj is defined only in the subspace generated by columnar 
states. The perturbation theory describing the effects of 
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small couplings oc J x acts in the subspace of columnar 
states^. We recall the expression of the effective Hamil- 
tonian obtained at leading order >&22 



H 



(0) 
col 



j=l...L* 



(2.8) 



Here the effective coupling constant J co \ describing the 
flip of a whole column is obtained at order L z in pertur- 
bation theory: 



J col = 2L z 2 L *jpJ* 



The coefficient 7P depends on boundary conditions 
(with or without a prime for OBC or PBC, respectively) 
and on the column length; it can be determined by con- 
sidering all processes flipping two neighboring columns i 
and (i + 1), by L z successive actions of perturbing terms 
J x o~fj<jfj +1 - Assuming PBC, the excitation energies of 
intermediate states during such L z -th order processes are 
integer multiples of the quantity 8| J z \ which appears in 
Eq. (|2.9[) . The counting of these processes, weighted by 
a factor depending on the excitation energies at each 
step of each process, is a combinatoric problem which, to 
our knowledge, does not have a general analytic solution; 
however for small L the exact values of 7l , or cquivalcntly 



8J Z 



(2.9) 



of c L = L2 L ~ 2 ^ y L ' , are easily computable. As examples 
we give here 74 = 5/4, 75 = 7/4 and 76 = 29/12. In the 
case of OBC, the number of processes flipping two neigh- 
boring columns of L sites is the same as for PBC, but 
the excitation energies at some intermediate steps may 
be lower than in the periodic case so that coefficients cl 
are be replaced by larger ones c' L , e.g. for L = 4 one has 
C4 = 32. With PBC one can even remark that > 1, by 
noticing that there are exactly L x 2 L ~ 2 column-flipping 
processes for which the excited energy at each step is 
minimal, i.e., 8\J Z \ (these are the processes where two 
successive actions of perturbing terms occur on bonds 
distant by 1 unit along the z axis). 

A scaling law for cl, or equivalently for the size- 
dependence of J C oi, was given in Ref. [H, indicating that 
the latter vanishes exponentially with increasing L z - 
in the compass model with \J X \ < \J Z \ this yields pre- 
cisely the 2 Lx -fold ground state degeneracy in the TL. 
The isotropic case J x = J z has a higher ground state 
degeneracy 2 Lm + 2 L = in the TL, which can somehow be 
deduced from similar arguments. 

Besides allowing to understand the unusual ground 
state properties of the QCM, the pseudospin formalism is 
also helpful in the framework of the CH model. There the 
low-energy spectrum of an adequate pseudospin Hamil- 
tonian consists of excitations characteristic of dominant 
compass interactions, the column-flips. These excitations 
which will be analysed in Sec. IVII correspond qualita- 
tively to flipping all spins in a column of a finite cluster 
[see Fig. Hid)]. They arc well defined when perturbing 
interactions favor a particular columnar pattern in the 



ground state, and we will see that this is typical for the 
CH model. 

As mentioned before, the compass model itself is char- 
acterized by a high level of frustration between a z o~ z and 
o~ x o~ x interactions, independent of the sign of the associ- 
ated amplitudes. From that perspective, the introduction 
of perturbing Heisenberg interactions seems to increase 
the degree of frustration in the model, e.g. in a case 
where they are of sign opposite to that of dominant com- 
pass interactions. The ordered patterns favored in this 
case, if ever, arc expected to differ from those selected for 
dominant Heisenberg interactions, i.e., \I\ J c - in the 
former case, a ground state minimimizing energy both of 
dominant compass interactions (on either vertical or hor- 
izontal bonds, depending on the sign of | J z \ — \ J x \) and of 
Heisenberg interactions (on other bonds) can be selected; 
while in the latter, a FM or Neel order is expected with 
an easy axis selected by compass couplings. Thus, be- 
sides the question of whether the exotic, semi-disordered 
ground states characteristic of the compass model can 
actually exist in presence of Heisenberg couplings with 
small amplitudes, one can focus in this model on the de- 
termination of the phase diagram, with multiple phase 
transitions between the more conventional FM or Neel 
phases, and more exotic phases. The characterization of 
these phases will be the object of the next chapter. 



III. SPIN STRUCTURE FACTORS 

In this Section we address the central question: What 
happens to the macroscopic 2 L -fold ground-state degen- 
eracy of the anisotropic QCM in the presence of Heisen- 
berg perturbations? We show that in the most general 
case, where compass coupling strengths | J x \ and | J z \ are 
non equal and where the Heisenberg coupling strength |/| 
is nonvanishing, the ground states are characterized by 
long-range spin order. This is evidenced by spin struc- 
ture factors S a (q), defined for each orthogonal spin com- 
ponent a G {x, y, z} as: 



(3.1) 



For given interaction parameters, a peak in S a (q) at a 
single momentum <fo and component a signals an order- 
ing of spins with a finite component along the a axis and 
a modulation period given by go- 

Let us exemplify this in a situation where all cou- 
plings are AF, and J z > J x . In Fig. 21 the large peak 
at M — (ir,ir) of S z (q) for / = 0.3J C indicates that the 
corresponding ground state is of the G z type, that is Neel- 
ordered with spins along the z axis. We encountered sim- 
ilar ordering features for I / J c varying from very small to 
very large values, keeping fixed the compass amplitudes 
such that < cj) < ir/4. In the I ^> J c limit, this can be 
interpreted as follows: small compass couplings break the 
SU (2) symmetry of Heisenberg interactions and select an 
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r x m y r 

FIG. 4: (Color online) Structure factors (a) S z (q) and (b) 
S y (q) (divided by N) obtained for the N — 36 cluster (sym- 
bols). The lines are guides to the eye. Parameters: cj> = 7r/10; 
the data correspond to four different values of I / '.J c : 0.3 (G z 
phase), +10 -6 , —0.3 (C' z phase) and —0.6 (F y phase). High 
symmetry points T = (0,0), A' = (tt,0), M = (tt, 7t) and 
Y — (0, tt) in the 2D Brillouin zone are defined in the inset. 

easy axis for Neel order; it is obvious here that this easy 
axis is z. In the limit /< J c , the selection of G z order is 
to be understood differently: we recall that compass cou- 
plings alone, on a L x x L z cluster, select a class of 2 L * 
columnar states, characterized by long-range Neel-typc 
correlations along columns but short-range correlations 
along rows. These states are separated from higher en- 
ergy levels by a gap which, in the Ising limit J x /J z 0, 
is ~ AJ Z . 

This semi-ordered nature of the GS is reflected in 
Fig. rjja) for I / J c = 10~ 6 by a structure factor spread 
over all momenta of the form q = (q Xl tt) - for J x = I = 0, 
one would indeed have S a (q) = L z 8 q _^. The effect of 
small J x compass couplings is mostly to reduce slightly 
the difference of S z (q) at {q x ,ir) and (q x ,q z ^ tt), re- 
spectively. In contrast, even very small Heisenberg cou- 
plings have a much stronger effect, seen in the example of 
Fig- HI they result in a strong enhancement of S z (tt,tt), 
compared to S z {q x ,ir) for q x ^ tt. The fact that this en- 
hancement is much stronger for I / J c = 10~ 6 , in this case, 
than for I / J c = 10~ 3 in the anisotropic case of Fig. [5jc), 
can be understood within the effective pseudospin model 
which we will describe in Sec. IVI A[ To explain this, we 
notice that the 2 ix columnar states include two Nccl-likc 
states, where spins on each x-orientcd bond connecting 
neighboring columns are antiferromagnetically arranged. 
These states are favored over the other 2 La: — 2 columnar 
states by AF Heisenberg couplings (even when I is in- 
finitesimally small) — namely by the o~ z o~ z components 
of these couplings on horizontal bonds. Transverse com- 
ponents of Heisenberg couplings have no matrix element 
between distinct columnar states. These favor an AF 
arrangement of the z component of spins on x-oricntcd 
nearest neigbor bonds, which at the global scale result in 



the G z order. We will provide in Sec. IVI Al an explana- 
tion, based on the formalism of pseudospins fj , for the 
fact that such an ordered phase can be selected in the TL, 
even for infinitesimal Heisenberg couplings. When these 
become larger, the G z order remains the most favorable, 
and is almost fluctuation-free for I < J z (evidenced in 
Fig. Ha), for I/J c = 0.3, by S Z (M)/N ~ 0.9 close to the 
maximal allowed value 1). 

A similar situation is also found for negative, moder- 
ate Heisenberg couplings I < 0, as shown in Fig. 0Jb) 
for I / J c = —0.3. Here, among columnar states, those fa- 
vored by Heisenberg couplings have FM correlations on 
x-oriented bonds. The selected ordered patterns define 
now a C z ordered phase (see cartoon in the correspond- 
ing region of the phase diagram for the AF CH model 
displayed in Sec. HVi which manifests itself by a peak at 
q = Y in the structure factor S z (q). 

With the same value of anisotropy parameter (f>, but 
large negative Heisenberg couplings (0 < J c -C —I), the 
selected order follows from a different mechanism. The 
dominant Heisenberg interactions tend to favor a FM 
phase, with possibly an easy axis due to the anisotropy 
induced by compass couplings. Actually, one sees in 
Fig. Bib) that S v (q) has a distinct maximum for / ~ 
— 0.6J C at q = r, while S z (q) (and S x (q), not shown) 
do not display such a peak. This indicates that the FM 
order develops, when increasing the Heisenberg coupling 
strength |/|, already for moderate values / = — 0.6J C , and 
has the easy axis y in spin space. Here, unlike in previous 
cases, a FM-ordered pattern cannot correspond to any of 
the columnar states favored by the compass interactions 
alone (since these have AF correlations on either x or z 
bonds), and the order along the y axis is a compromise. 

The appearance of the i^-phase can be understood as 
follows: The AF compass couplings frustrate the dom- 
inant FM Heisenberg ones; but since they act only on 
two spin components, the system avoids this frustration 
by rotating spins from the xz plane to the y axis where 
they fully profit from Heisenberg interactions. Although 
this mechanism is based on a classical picture, it explains 
the behavior observed over a wide range of I / J z includ- 
ing the example shown / = — 0.6J 2 where frustrating 
compass terms are of amplitudes comparable to Heisen- 
berg ones. Several other ordered phases, displayed in 
the phase diagrams shown in Sec. IIV[ are found when 
at least one of the coupling constants J x , J z is negative. 
For small \I\/J C , they always result from the selection of 
a pair of compass states by small Heisenberg couplings, 
because they couple components along the easy axis [z 
or x) of spins neighboring on x or z bonds, respectively. 

We turn now to the case where compass couplings are 
isotropic, i.e., J z = J x . There, in absence of Heisen- 
berg couplings the low-energy states of a L x x L z sys- 
tem consist not only of 2 L * columnar states, but also of 
2 Lz row-type states The latter ones have spins oriented 
along x and long-range correlated along rows, but not 
along columns. As shown in Ref. [H this situation corre- 
sponds to a first order transition point between two dis- 
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FIG. 5: (Color online) Structure factors (a-c) S z (q) and (d- 
f) S y (q) obtained for the isotropic compass model J x — Jz 
(4> = 7r/4) and increasing values of the Heisenberg interaction 
I for clusters of different size: (a,d) N = 16, (b,e) N = 24, 
and (c,f) N = 36 clusters. The high symmetry points V, X, 
Y and M are as in Fig. [3] The compass quasi- ID order is 
characterized by constant values of S z (q) along the M — Y 
line. The onset of AF order in the TL is evidenced by the 
increasing maximum of S z (q) (it, n) which is visible already 
for very small values of I > 0. 



tinct phases characterized by either column- (J z > J x ) 



or row- (J z < J x ) -type ground states in the TL. Yet, 
as in the previously discussed anisotropic case, in the 
isotropic one J x = J z small Heisenberg couplings select 
among these states only a small number, here four. The 
selected states are here Neel states: two of them have 
spins along z, while the two others with spins along x 
are selected within the class of row-type states. These 
four states are a priori characteristic of a Z2 x Z2 or- 
dered phase, breaking spontaneously not only transla- 
tion symmetry, but also the symmetry {cp, erH, <t5} — > 

}, where R is the reflection w.r.t. 
the z = x diagonal in real space. 

The Z2 x Z2 ordered phase is characterized by two 
peaks in structure factors S x (q) and S z (q), both at 
q = M. Structure factors S z (q) are shown in Figs.^a)- 
[5Jc) for 3 different clusters of N = 16, 24, 36 sites and 
several values of I / J c ranging from to 5. In contrast 
with the anisotropic case where such a peak can have a 
value close to the maximum allowed (N, number of sites), 
here peak amplitudes are limited by sum rules to N/2 on 
isotropic clusters with L x = L z [Note that if the cluster 
is anisotropic the peak amplitudes can slightly exceed 
the value N/2, as in Fig. E^b) for (L X ,L Z ) = (4,6) and 
I/J c = 0.1]. A more striking feature is that these peaks 
grow very fast with small /: peak amplitudes exceeding 
75% of the maximal value are attained for / = 0.01J C 
on the largest clusters. In a situation with a slight com- 
pass anisotropy \J X — J z \ <C J z one can show (see Ref. 
IM ) that the order develops as soon as I ^> J co \ with 
Jcoi vanishing exponentially with increasing L z ; this ar- 
gument extends to the isotropic case (see also Sec I VI Al 
for details). In contrast, we also show structure factors 
S v (q) in Figs. EJd)-[SJf). A peak is also observed when 
Heisenberg and compass couplings are of similar values, 
but the peak amplitude is almost independent of system 
size (consider e.g. the I = J c case). For these reasons one 
can conclude that the Neel order, with spin directions x 
and z equally favored over the y direction, is selected in 
the TL for arbitrarily small Heisenberg couplings, in the 
whole range of values // J c > in the isotropic AF case. 
We will see in the next Section that this order is unstable 
even to infinitesimal variations of J z — J x : depending on 
the sign of this quantity, either the Neel patterns with 
spin along z or those with spins along x are favored, and 
one recovers the G z or G x phases. 



IV. PHASE DIAGRAM 

The CH model reveals a large variety of ordered ground 
states as function of the interaction parameters J x /J z 
and 1 1 J c . Some of these phases were described in the 
previous Section. The determination of the ground state 
phase diagram, and the characterization of QPTs (as, for 
instance, the G z — G x transition discussed above) will be 
the object of the present Section. We will first give a 
classification of the possible phases expected in the clas- 
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sical limit of the model; then we will determine the phase 
diagram, first restricting ourselves to the AF CH model 
(case J z > 0), before addressing properties specific to the 
FM CH model (case J z < 0). 



A. Ordered phases of the CH-model 

To analyze the phase diagram of the compass- 
Heisenberg model, we will consider first the classical (or 
large S) limit, where one regards spins as vectors living 
on a unit sphere. This is of prime interest, since we will 
see that all ordered phases of the model are found in this 
limit. To draw a tentative classical phase diagram one 
needs to compare the ground state energies Eq associated 
to these different phases. In Table U we present a list the 
candidate phases $ Q . For each phase the index a G x, y, z 
denotes the easy axis or spin direction favored, while the 
capital letter in <£>„ indicates the type of spatial structure 
or correlation pattern, i.e., G for Neel-type AF phase, F 
for FM phase, and C for columnar or C-type AF order, 
i.e., with nearest neighbor spin correlations AF for one 
bond direction and FM for the other. By convention, the 
presence of an index in $ for C-type phases indicates 
that nearest neighbor correlations are AF on bonds where 
compass interactions couple spin components along the 
easy axis. Furthermore an ordered state $ is determined 
by the momentum q for which S a (q)/N stays finite in the 
TL in this phase. Eventually we give the ground state 
energy (per site) Eo($ a ) of each phase in the classical 
limit. 

For instance, the F z phase is the FM phase with spins 
along the z axis - this order implies that the structure 
factor S z (q) develops a peak at q = L of amplitude fi- 
nite in the TL. The classical energy per site is the sum 
of a contribution of compass couplings on z bonds, and 
of Hciscnberg couplings on all bonds. The C x phase in- 
stead is an ordered phase with spins along x, so that the 
ordering wave vector (at which S x (q)/N is finite in the 
TL) is Y = (0, tt), i.e., with FM spin correlations along x 
bonds and AF along z bonds. In the classical limit only 
compass couplings contribute to its energy per site J x , 
since the contributions of Hciscnberg phase on x bonds 
and on z bonds cancel each other. 



TABLE I: Classification of ordered phases of the CH model. 
For each phase $, the easy spin axis a, the ordering wave vec- 
tor q and the classical energy per site Eo{&) are given. Here 
$ indicates either a FM (F), Neel-type AF (G) or columnar 
(C) phase. 
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It is clear that each of these phases is stabilized, at 
least in the classical version of the model, when the cou- 
pling amplitude(s) entering Eq(§) is/are much larger in 
absolute value than other amplitude (s) - for phases F x / Z 
and G x i z the condition is that both amplitudes involved 
in Eo(<&) have equal sign and that their sum is much 
larger than the amplitude not involved. The determina- 
tion of the domains of stability of these phases, first in 
the classical limit and then in the S — 1/2 model, will be 
described hereafter, focusing first on the AF case J z > 0. 

B. Phase diagram: antiferromagnetic case J~ > 0. 

1. Classical approach and symmetry relations 

Here we consider the case J z > 0, where most of the 
phases listed in Table U — actually all but F z , C z and G y 
can be stabilized depending on the values of interaction 
parameters I / J z and J x / Jz- By using the classical en- 
ergies in this Table and determining, for fixed I / J z and 
Jx/ Jz, which of these energies is the lowest one (with the 
assumption, which will appear in the following as justi- 
fied, that no other phase is stabilized in a finite volume 
of the phase space determined by these two parameters) , 
we find the classical phase diagram represented in Fig. |6l 
with phase transitions corresponding to continuous lines. 
The classical phase boundaries are straight lines in the 
present parametrization, because the classical energies 
depend linearly on the various coupling amplitudes. It 
is not surprising to see that for all interactions being AF 
(Jx,Jz,I > 0) Neel order is always favored, with a G z 
or G x phase depending on the sign of J z — J x . More 
interesting is the extent of the C' z phase for moderate, 
negative /. Although in this phase, due to the z orien- 
tation of spins, the J x compass couplings are frustrated, 
their sign matters for the stability of this phase: its ex- 
tent in the phase diagram is smaller for J x < than for 
J x > 0, since in the former case it competes with a F x 
phase stabilized by J x couplings, whereas in the latter 
case it competes with a F y phase where J x couplings are 
inactive (which explains why the C' z — F y transition line 
is independent of J x /J z ). 

Among these phase transitions, several ones have a po- 
sition in the phase diagram of Fig. [6] which can be deter- 
mined from symmetry considerations, and is thus insen- 
sitive to quantum fluctuations. The simplest example is 
that of the G x — G z transition: intuitively, one can guess 
that it can occur only for J z = J x (and I > 0), but one 
can also notice that the transformation defined by: 

{T/,T s \T/) = {a%*l-CT% (4.1) 

where s = R(r) and R is a spatial rotation of ir/2 around 
a reference site, allows to rewrite the CH Hamiltonian as 
follows: 

H = J x ^^ TijT^ij + J z ^ T^j^j+l 
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FIG. 6: (Color online) Phase diagram of the Compass- 
Heisenberg model in the ( J m , J)-plane for fixed AF interaction 
J z = 1. Long-range order is stabilized by any finite I. Square 
(J x = J z ) and diamond (J x = — J z ) indicate multi-critical 
points, where in each case four ordered phases meet. The 
spin order of the different phases {G z , G X ,C' Z , C X ,C' X , F x , F y } 
is depicted in a corresponding inset, and the subscript a = 
x, y, z indicates the type of symmetry breaking in spin space. 
The QPTs between F x and C' z phases, and between C x and 
G z phases (solid lines) are modified by quantum corrections 
w.r.t. the dashed straight lines which show the corresponding 
classical transitions. 

+ ^^•■(^■+1+^). (4.2) 

In other words, this transformation maps the domain 
Jz < Jx of the phase diagram onto the domain J z > J x 
and vice- versa, and if a point with J z < J x is in the G z 
phase it implies that its image by this transformation is 
in the G x phase. Only at the transition line J x = J z the 
CH Hamiltonian is invariant when this transformation is 
applied which means that the G x — G z transition has to 
occur there (unless another intermediate phase is stabi- 
lized, which is not expected). For I < < J X ,J Z , the 
same transformation is a bijcction between each point of 
the C' z phase, with given I/J z and <j) = tan _1 ( J x j J z ), 
and a point in the C' x phase, with the same value of I / J z 
and anisotropy parameter 7r/2 — <fr. This implies that the 
transition line between these two phases is fixed to the 
line J x = J z as in the classical limit, under the condition 
that the |/| is too small to stabilize the FM F y phase 
which would be favorable otherwise. One can notice — 
here at the classical level but this feature is actually con- 
served in the quantum model — that the isotropic point 
of the compass model, with J x = J z > and I = 0, is 
unstable to even infinitesimal variations of either J x — J z 
or I; depending on the sign of both quantities four differ- 
ent phases can be selected, so that this point can be seen 
as a quadricritical point in the context of the CH model 
(we do not mean by this that correlations are algebraic 
there, but simply that four phases meet at this point). 

Another transition characterized by an additional sym- 
metry is the one between the C' z and the F y phases, ob- 
tained by varying the Heisenberg amplitude and keeping 



4> € (0;7r/4) fixed, and stabilized for < —I <C J c and 
for —I 3> J c , respectively. In this case one can make use 
of a transformation defined by: 

{^,T{j,r;;;} = {(-l)V^,^,(-l)V^}. (4.3) 

Reexpressing all couplings in function of r' operators, one 
obtains the following Hamiltonian: 

H = H x [t" x ] 

- (i + Jxtfj-Ti+ij + . (4.4) 

Here H x [t x ] is a function of only x-components of t" 
spins, necessarily invariant under further rotations of 
spins along the x axis (not combined, here, with any 
spatial symmetry but the identity) . From the expression 
of H - H x [t" x ] in Eq. one sees that for I = -J z /2 

such rotations leave H invariant: there, an extra U(l) 
symmetry appears, so that all states with fully polarized 
t" spins in the yz plane are degenerate ground states. 
These contain the two degenerate ground states of the 
F y phase, as well as those of the C' z phase (which is also 
FM in terms of t" spins): necessarily, the transition be- 
tween both phases has to occur there. Notice that, unlike 
for previously discussed G x — G z and C' x — C' z transitions, 
here no mapping from the F y to the C' z phase is allowed 
away from the transition line. 

In contrast, another phase transition occurring in the 
classical phase diagram found in Fig. [6] for I < and 
— J z < J x < 0, namely the one between C' z and F x 
phases, is not fixed by any symmetry relation. Not only 
one cannot find a mapping between both phases as in the 
G z — G x case, but also if one looks for a transformation 
of the type given by Eq. (|4.3p . there is no way to iden- 
tify an extra symmetry (t/(l) or other) at the classical 
transition line / = — (J z + J x )/2. In consequence, the 
corresponding transition line can be shifted by quantum 
fluctuations — which of the two phases is stabilized by 
those fluctuations at the I = — (J z + J x )/2 line is one of 
the questions addressed in the next paragraph. 

2. Phase transitions in the quantum CH model 

We turn now to the phase diagram of the S — 1/2 CH 
model and investigate those aspects which show up as a 
result of quantum fluctuations. In Sec. IIV Al the mag- 
netically ordered phases were selected either for |/| <C J c 
or |/| 3> J c , and subsequently their respective stability 
in the classical phase diagram was discussed. There we 
made the implicit assumption that no other phase occurs 
in an intermediate range of /; we will now see that this 
assumption is justified even in the quantum model. 

The first case to be addressed is the transition between 
the C' z and F y phases, occurring by increasing \ I\/ J z from 
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FIG. 7: (Color online) Evolution of order parameters 
S v (0,0)/N and 5 z (0,7r)/iV across the C' z - F y transition, by 
varying I / J z ; anisotropic compass interactions have the same 
parameter <f> = 7r/10 as in Fig. [4] The cluster N = 24 has 
dimensions = 4, = 6. The dashed-dotted lines cor- 
respond to perturbative estimations of order parameters — 
given for the C' z phase by Eq. (|4.5|) . The transition coincides 
here with the classical value I c = Jz/2. 



to infinity, with I < < J z and fixed J x /J z € (0;1). 
The classical approach (see Sec. IIVB 1[) predicts a tran- 
sition between those phases at I = —J z /2. In the case of 
the quantum model we study in Fig. [7] the evolution of 
the spin structure factors S Z (Y) and S V (T) correspond- 
ing to the C' z and F y phases as function of I/J z . The 
data in Fig. [7] calculated at fixed <fi — tt/10 indicates that 
no intermediate phase is stabilized in a finite range of I 
between the C' z and F y phases, since on both sides of the 
classical transition point I c = —J z /2 (I c ~ — 0.4755J C ), 
cither S Z {Y)/N or S V (T)/N takes large values. This is 
clear evidence for long-range magnetic order of either C' z 
or Fy-type, respectively. The transition is clearly de- 
tectable already on clusters of moderate size, by the sharp 
evolution in its vicinity of structure factors as function 
of 1/ J z : the maximal slopes are found exactly at I = I c . 
The size scalings of both order parameters (i.e., of the re- 
lated structure factors divided by N) arc shown in Fig. [8] 
and provide evidence of this transition on a more quan- 
titative level, with for each order parameter exhibiting a 
clear change of behavior at / = I c : for |/| < |J C | the scal- 
ings indicate that S Z (Y)/N and S V (T)/N take in the TL 
respectively a finite value or 0, while for |/| > \I C \ it is the 
contrary — eventually in the particular case I = I c both 
order parameters scale down to a common finite value, 
confirming as well that this transition point can also be 
seen as an intermediate phase where a U(l) symmetry 
is spontaneously broken and the ground state is an XY- 
type ferromagnet in terms of r" spins defined in Eq. (|4.4p . 
The corresponding order parameter, (S V (T) + S Z (Y))/N, 
would be equal to 1 in the TL if terms H x [t x } were absent 
from Eq. (|4.4j) : in their presence, for <f> = tt/10, Fig. [^in- 
dicates that this order parameter attains ~ 0.7(1) in the 



FIG. 8: (Color online) Size-scaling of order parameters of (a) 
C' z and (b) F y phases for fixed anisotropy parameter cj> — rr/10 
and for values of I / J c at and close to the value / = —0.4755 J c 
of the transition between both phases. 

TL. Yet this XY-type order is sensitive to infinitesimal 
variations of /, which lower the symmetry of the Hamil- 
tonian from U(l) to Z2 and select an easy axis y or z for 
the ordered phase. 

Far away enough from this phase transition, quan- 
tum corrections to the values of order parameters (which 
would show, in the classical limit, a discontinuity at I c 
with a jump/drop from 1 to or from to 1) arc relatively 
well estimated using second-order perturbation theory^ 
There, for a given phase the unperturbed Hamiltonian 
consists of components along the easy axis (e.g. z in the 
C' z phase) of both compass and Hcisenberg couplings, 
while transverse components of these couplings are re- 
garded as perturbations. One can illustrate this in the 
case of the C' z phase: the effects of quantum fluctuations 
is to reduce (in absolute value) the correlation between a z 
components of spins, and thus the C' z order parameter. 
For sites r= and s = situated at distance 

d > 1 from each other, one obtains: 

x { 1 + W^ny + (J^i?} '• (4 ' 5) 

The prefactor (— 1) M cancels out with the phase fac- 
tor in Eq. (|3.1|l giving the order parameter S Z (Y)/N, 
which in the TL is equal to the absolute value of the ex- 
pression in Eq. (|4.5p (corrections on finite systems result, 
within this approach, from the contributions of nearest 
neighbor and on-site spin correlations: only at short dis- 
tances d < 1 the correlations differ from the long-distance 
limit, i.e., entanglement between further neighbor spins 
is neglected at this order in perturbation). The pertur- 
bative estimate >S z (Y;pert) of the C' z structure factor is 
shown in Fig. [7] for cj) = 7r/10, and gives good agree- 
ment with the numerics, away from the transition to the 



11 



s°Vn 




FIG. 9: (Color online) Variation of S x (0,0) and S z (0,tt) 
across the C'~ — F x transition as function of I/Jz at fixed 
<j> — — 37r/20, and for system sizes N = 16, 4 x 6, 32. Here the 
quantum phase transition is shifted by quantum fluctuations 
from I c = -0.245 J 2 to I c = -0.204 



F y phase (here, for |/| < 0.3). Similar estimates can 
be done in other ordered phases, like S y (T; pert), also 
shown on Fig. [JJ in the F y phase: here the unperturbed 
Hamiltonian consists of couplings cEcH in T-L. Concern- 
ing absolute values, the agreement with numerics is less 
accurate than in the C' z phase but the dependence on 
1 1 ' J z is correctly reproduced by the perturbative result. 

The spin structure factors are even more useful to 
study transitions not characterized by additional sym- 
metries, such that the phase boundaries can be modi- 
fied by quantum fluctuations. As an example we focus 
on the C' z — F x transition: the relevant order param- 
eters are S Z (Y)/N and S X (T)/N. We show in Fig. 1 
their evolution as function of I/J z , again for fixed J x j J z . 
For each cluster size, the two curves have maximal slope 
at the same value of I/J z , which can thus be consid- 
ered as a finite-size transition point. But in contrast to 
the C' z — F y case, here this transition point is cluster- 
dependent and distinct from the classical transition point 
(7° ~ -0.245J 2 for = -3tt/20). The dependence on N 
of this finite-size transition point is rather weak, and this 
allows us to locate approximately the transition point in 
the TL — in the case of parameters of Fig. it occurs 
at I c /J z ~ -0.21(1). 

The deviation of this transition point from the classi- 
cal value I® can be well estimated by evaluating energies 
of both phases using second-order perturbation theory. 
This approach gives the following estimates for the ener- 
gies per sites of different phases: 



E(F X ) = E (F X ) + 
E(C' Z ) = E (C' Z ) - 



J 2 



8J X + nr 

8J Z + 41 J z - 



(4.6) 
(4.7) 
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FIG. 10: (Color online) Evolution of 1 — /, where / denotes 
the fidelity, as function of 1/ J z for parameters (a) (j> — —3n/20 
and (b) cf> = 7r/10 and for cluster sizes N = 16 and iV = 24 (in 
both cases L z =4). Smooth peaks reflect phase transitions 
between the F x or F y phase and the C' z phase that are con- 
tinuous for finite systems. The sharp peaks at I = indicate 
the C' z — G z transition. 



(2) 

case 4> = — 37r/20, this value is Ic — — 0.2040 J c , that is 
very close to I c estimated from order parameters of both 
phases. More generally, for variable J x /J z € (— 1;0), the 
transition line between C' z and F x phases as estimated 
from Eqs. (|4.6[) and (|4.T[) is shown on the phase dia- 
gram, see Fig. [51 One sees there that the deviation from 
the classical transition is always in the same manner, i.e., 
at the classical line I = I® = — (J z + J x )/2 quantum fluc- 
tuations favor the F x phase at the expense of the C' z 
phase. This may appear surprising since — at least in 
unfrustrated Heisenberg models — a FM phase allows for 
fluctuation-free ground states, while orders with some AF 
bonds, like the C' z phase, are typically accompanied by 
quantum fluctuations. The case of the C' z — F x transi- 
tion is different: first, both phases are characterized by 
easy axes, distinct from each other; and second, contribu- 
tions from different bonds to quantum fluctuations have 
to be considered separately. The contribution of z-bonds 
to quantum fluctuations, thanks to the large amplitude 
Jz > \Jx\, \I\ in the vicinity of the classical transition, 
removes the degeneracy and stabilizes the F x phase with 
respect to the C z one. 

Complementary to the detection from structure factors 
of the C' z -F y and C' z -F x phase transitions, one can also 
analyze the behavior as function of I / J z of the fidelity^ 
defined as: 



/(J) = K*o(/ + W)|*o(J-M)>|, 



(4.8) 



Within this approach, the transition point is given by 
the value of I for which E(F X ) = E(C' Z ): still in the 



where both GSs are computed for values of Heisenberg 
amplitudes differing by ±51 (51 = 5 x 10~ 3 J c ) from the 
nominal value I. In Fig.QUlwe plot the quantity ln(l — /) 
as a function of I/J z for two clusters (with N = 16 or 
N = 24 sites), and for the two values of <j> correspond- 
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FIG. 11: (Color online) (a) Lowest excitation energies AE 
across the F x — C' z and C' z — G~ transitions as function of 
1 1 J z at cj> = — 37r/10 and for different cluster sizes; (b) Ground 
state energy Eo reflecting the avoided crossing of the F x — C' z 
transition at I/J z = —0.21(1). 



ing to Figs. [7J and [5J respectively. For a given cluster 
size and value of </>, peaks are observed on the I / J z axis 
at positions coinciding with the maximal slopes of order 
parameters in Figs. [7J and These peaks are thus good 
indicators of phase transitions, here between the C' z and 
o ther phases. Note that the peaks at / = (transition 
between C' z and G z phases on the compass line of the 
phase diagram, either for J x > or J x < 0) are much 
higher and thinner than those at transitions between the 
C' z phase and either the F y or the F x phase, for > 
and 4> < 0, respectively. Indeed, the qualitative change 
in the GS occurs continuously at the C' z — G z transi- 
tion, but in a very narrow range of I / J z (estimated in 
Sec. I VI Ap . resulting in a sharp peak centered on / = 0. 
In contrast, at the transitions between the C' z and FM 
phases the peaks in ln(l — /), centered around I c (up 
to small deviations resulting from finite size), are much 
more smooth, characteristic of a continuous transition. 
The latter behavior may be, a priori, an artefact which 
follows from finite size, while we have indications due to 
sharpening of peaks with increasing system size that the 
transition may become first order in the TL. 

Eventually, the phase transitions in the CH model can 
also be addressed by considering the low-energy spec- 
trum, which we do once again on the example of the 
C' z — F x transition for fixed <fi = — 37r/20. The depen- 
dence of the ground state energy on I / J z , shown in Fig. 
[TTT down) . is consistent with that of the fidelity: Eo/N for 
a given cluster varies smoothly in the vicinity of the tran- 
sition, but from the comparison between different cluster 
sizes it appears that a cusp develops in the TL at J = I c , 
which supports the picture of a first order transition (in 
this limit) responsible for the peak of ln(l — /). The 
nature of the transition may also be examined by con- 
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FIG. 12: (Color online) Lowest excitation energies AE for 
N = 16 across the transitions between the F y and C' z phases 
(for three distinct values of J m /J z between and 1, corre- 
sponding to (j> — 7r/10, 7r/5, and 7r/4). 



sidcring the lowest excitation energies, see Fig. [TTt up). 
On each side of the transition, the lowest excited state is 
found in a representation indicative of the symmetry of 
the phase: 

(a) In the F x phase, it is found, as the ground state, 
with k = T; but both states have opposite parity of 
^_,er£/2 (the latter quantity being conserved in the 
model). The energy splitting between both states 
increases when approaching the transition point, 
and for fixed I / J z decreases (exponentially, as far 
as one can tell) to zero with increasing linear size. 
This means that in the TL both states become de- 
generate ground states; among linear combinations 
of them one finds states fully-polarized either along 
+x or —x in spin space, and which break sponta- 
neously the global symmetry under cr5 — > — cr£. 

(b) In the C' z phase, the second lowest state, also with 
an excitation energy decreasing rapidly to zero with 
increasing size, has identical parity of h J^? crS as 
the GS; but a distinct momentum Y = (0, ir). Here 
these states, degenerate in the TL, are characteris- 
tic of the C^-type order, with spontaneously broken 
symmetry of translation by one lattice unit along z 
bonds — or, in terms of symmetries in spin space, 
spontaneous breaking of the global symmetry un- 
der a- — ► — takes place. 

At equal size, the positions of the crossings seen on 
Fig. [TT] for the C' z - F x transition, and Fig. [TJ for the 
C' z — F y transition, match well the positions of maximal 
slopes in Figs. M an d UJ respectively, and at the cross- 
ing their common excitation energy seems to decrease to 
zero towards the TL. Such level crossings on finite sys- 
tems are thus good indicators of the corresponding phase 
transitions. 
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Considering the various features of phase transitions 
described above, we can distinguish several types of tran- 
sitions. Those occurring at I ^ 0, such as between C' z 
and either F y or F x phases, require a particular atten- 
tion. Although no level crossing is observed in the ground 
state on finite systems at these transition, several fea- 
tures indicate that they might have a first-order charac- 
ter in the TL: (i) the ground state energy, as function of 
the interaction parameter driving the transition, seems 
to develop in the TL a cusp characteristic of a first or- 
der transition^ (ii) the ordered phases on each side have 
distinct Z2 symmetry groups (or distinct spontaneously 
broken symmetries); (iii) accordingly, a crossing occurs at 
the transition between the two lowest excitations, found 
in different symmetry representations, and each of these 
excitations becomes in the TL, in one of the two phases, 
i.e., one of the two degenerate ground states of this phase; 
(iv) tentative scalings of order parameters suggest that 
in the TL they jump at the transition between zero and 
a finite value. 

However, these indications are no evidence for a first 
order character as scalings may be biased by the small 
system sizes available. More importantly, the fact that 
the two competing phases have distinct symmetries does 
not prohibit a continuous transition, although beyond 
the Landau-Ginzburg paradigm, between these phases^ 
Here, the vanishing of the lowest excitation energy can 
also signal that the system becomes becomes gapless at 
the transition. This is clear in the C' z — F y transition 
along the / = — J z /2 line, where the U(l) symmetry 
is spontaneously broken: the finite values in the TL of 
S Z (Y) and S V (T), which in the rotated basis of f spins 
are here two components of the order parameter for an 
XY-type ferromagnet. Schematically, by varying I/J z 
the ordered moment can be rotated continuously from z 
to y at the transition point, in contrast to typical first 
ordered transitions where hysteresis phenomena usually 
occur. More generally, we will see in Sec. [V] that each 
phase transition away from the 1 = line is characterized 
by the vanishing of the anisotropy gap to spin waves, 
which is finite in the Z2 ordered phases on each side of 
the transition. In consequence, the hypothesis that these 
transitions are continuous not only on finite systems but 
also in the TL is likely as well. We also remark that, 
whereas some of these transitions (as the C' z — F y one) 
are particular, with an an additional symmetry at the 
transition line, the same features occur at transitions not 
characterized by such additional symmetry. 

Eventually we comment on the 1 = line of the phase 
diagram, which can be seen as a transition line between 
distinct ordered phases (e.g. between the C' z and G z 
phases by increasing / from negative to positive values). 
There, spin waves are gapped (except at isotropic points 
where \J Z \ = \J X \)', these transitions are not character- 
ized by the softening of spin waves, but of column-flips 
introduced in Sec. ILTland which are gapless in the TL for 
1 = 0. This transition line, where one recovers the com- 
pass model, is characterized by the non-local invariants 
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FIG. 13: (Color online) Phase diagram of the compass- 
Heisenberg model in the (J x /J z , i/J z )-plane for fixed FM 
coupling J z = —1. Similar to Fig. [5] square (J x = J z ) and 
diamond (J x = —J z ) indicate multi-critical points, and the 
spin order of the different phases is depicted in insets. Phase 
labels follow the same convention as in Fig. [U The QPTs 
between F z and C' x phases, and between C z and G x phases 
(solid lines) are modified by quantum corrections w.r.t. clas- 
sical transitions (dashed lines). 

of Eq. (|2.6p . and the evolution between two distinct Z2 
ordered phases through this line can hardly be classified 
as an usual first- or second-order transition. 



C. Phase diagram in the ferromagnetic case J z < 

In the previous paragraphs we restricted to the case 
J z > and described the corresponding phase diagram 
and phase transitions, in function of two interaction pa- 
rameters: J x / J z and I/J z . Here we address the comple- 
mentary case with FM couplings on z bonds, i.e., J z < 0. 
The phase diagram of the ferromagnetic CH model is 
shown in Fig. 1131 It has many similarities with that of 
the AF CH model in Fig. [BJ There is an obvious differ- 
ence in the nature of the various phases, with e.g. for 
Jx/Jz G [0; 1] and I / J z > a F z phase replacing the G z 
one of the J z > case. Another qualitative difference be- 
tween both cases is that — here we assume \J Z \ > \ J x \ - 
the column-ordered states which are favored by dominant 
compass interactions allow for significantly less quantum 
fluctuations in the present J z < case than in the J z > 
one. Concretely, this comes from the fact that Heisen- 
berg terms on vertical bonds are inactive on columns with 
all spins aligned; this difference matters for the structure 
of the low-energy spectrum. Eventually, another moti- 
vation to study the J z < case follows from a possible 
qualitative description of arrays of NV centers coupled 
by quasi-short-ranged dipolar interactions [see Eq. (jl.ll) ]: 
the situation most likely captured within the CH model 
is with FM compass couplings J XtZ < 0, accompanied by 
AF Hciscnbcrg couplings. 

We address here the main ground state properties of 
the FM CH model by considering first the classical limit. 
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There, its phase diagram actually has the same topol- 
ogy as in the AF case: transition lines have the same 
positions and only the nature of the long-range order 
in each individual phase, stable in a particular range of 
{J x / Jz,I / Jz\ parameters, is determined by the sign of 
J z . This is because in this limit and due to the ab- 
sence of interactions between spins of the same sublatticc 
(in terms of bipartite or Neel sublatticcs), a transforma- 
tion inverting the sign of all couplings and simultaneously 
changing spins a? to — ap on one sublattice leaves the en- 
ergy in the classical limit unchanged — thus the ground 
states of both phases under consideration are related by 
a spin reversal on one sublattice. 

Coming back to the quantum model: here the shapes of 
the phase diagrams of the FM and AF model differ from 
each other in the same regions where they differ from 
their respective classical counterparts, since the transi- 
tion lines which are not fixed for symmetry reasons are 
differently affected by quantum fluctuations in the J z < 
case than in the J z > case (Fig. [5]). The case to com- 
pare to the previously discussed C' z — F x transition is here 
the C z — G x transition, which classically occurs on the 
line I = — (J z + J x )/2 for < J x < —J z . Here as well, 
one can estimate in second order perturbation theory the 
energies per site of both phases: 



Holstein-Primakoff bosons {at}, defined as follows, 



Sp = S- ata f = 5 - n P , Si = V2S]^ 1 - ^ a?, 

(5.1) 

with the usual spin-1/2 operators S? = \op. Due to the 
lack of SU(2)-invariance of the model, after linearization 
the spin- wave Hamiltonian contains not only ata^type 
but also a?ag-type terms, that do not conserve the num- 
ber of bosons: 

#LSW = 45 ^2 { [ la ?( a ?+?. + a r+e, ) + H.c] 

r 

-(4J + 2J z )np+ -J x (ataft + g x + apap +Sx + Jf.c.)}{5.2) 

with e x = (0,1) and e z = (1,0). Therefore, and un- 
like the nearest neighbor FM Hcisenberg model, the spin- 
wave dispersion does not depend linearly on coupling am- 
plitudes { J x , J z , I}, but as in the AF Heisenberg case has 
a square-root form: 



E(C Z ) = J z 



{2I + J X ) 2 

81 JJ -4:1' 



(4.9) 



From this one finds that on the line of the classical phase 
transition the energy of the C z phase is lower due to 
quantum fluctuations than that of the G x phase. This 
implies that, when taking quantum fluctuations into ac- 
count, the C z — G x transition line in Fig. 1131 must have 
a curvature opposite to the C' z — F x line in Fig. |6] Here 
again, the perturbative estimation of transition points, 
which yields for example I c /J z = —0.295 at J x /J z = 
—0.5, matches well with the numerical estimates from 
data of structure factors (not shown) of finite clusters. 



V. SPIN WAVE EXCITATIONS 

In the previous Sections we have mostly focused on 
ground state properties of the CH model, and considered 
lowest excitation energies merely as a tool to character- 
ize the symmetry of ordered phases and to locate phase 
transitions. Here we provide a description of the lowest 
excitations characteristic of the various ordered phases in 
the TL, these excitations being spin waves; for this we 
will use linear spin-wave (LSW) theory and see that this 
describes efficiently the lowest single-magnon branches. 

We begin the analysis of spin waves with the case of a 
FM phase, namely the F z phase corresponding to J z , I < 
and \J X \ < \J Z \. The classical ground state with all 
spins pointing along +z corresponds to the vacuum of 



co Fz (k) = 45^(21 J z + 2I\ + J kx ) - J 2 kxl (5.3) 
{Jkx-,1%} = {JxCOsk x ,2I(cosk x + cosk z )}. (5.4) 

The dispersion Eq. (|5.3j) is shown in Fig. [Mt a) on a 
closed path r->l->Jl/->y^rin the first Bril- 
louin zone for the coupling constants: I = —0.2J C and 
(f) = 7r + 37r/20. The LSW approximation describes well 
the lowest excitation energies obtained by exact diago- 
nalization (shown in the same figure for several periodic 
clusters). In diagonalization spin waves can be identified 
by the parity of J^f^p * na ^ ^ s opposite to that of the GS. 
One notices that for momenta such that ui(k) exceeds a 
certain critical value, spin-wave energies obtained by ex- 
act diagonalization tend to differ from the LSW results. 
Furthermore, this critical value seems to increase with 
system size. These features are actually related to the 
presence, in the low-energy spectrum of finite clusters, of 
the column-flip excitations which will be analyzed in the 
next Section. When the energy oc L Z \I\ of such excita- 
tions coincides with that of spin waves, their interaction 
leads to deviations from the LSW dispersion. Note that 
in general the spin waves obtained in the LSW theory 
are gapped, except in two limits: J x j I , J z /I — > (one 
recovers here the dispersion of the n.n. Heisenberg fer- 
romagnet), and for J x = J z . The softening, occurring at 

k = r, is associated to the transition to the F x phase. 

We turn now to the description of spin waves in AF 
phases. Here, the vacuum of Holstein-Primakoff bosons 
has to be defined differently; one can apply a transforma- 
tion similar to those seen in Sec lIVPTTl that is, inverting 
two components of spins on one sublattice, in order to ob- 
tain (after this transformation) a Hamiltonian with FM 
classical ground states. A concrete example is given here 
with the C z phase found in the phase diagram Fig. [131 
Here such a transformation consists of inverting y and 



15 




- LSW 



X M x y r 
k 




O N=32 
- LSW 



°T X M x y~ r 

k 




O N=32 



°r x Hvi x y r 

k 

FIG. 14: (Color online) Spin wave dispersions obtained by 
LSW theory (lines) and by exact diagonalization (symbols) 
of various clusters with size N < 32 for (f> = 237r/20 and for 
different Heisenberg coupling constants: (a) / = — 0.2J C ; (b) 
I — 0.2J C ; and (c) I = 0.6J C . These different parameters 
correspond to the F z , C z and G y phases, respectively. The 
abscisse Xk follows a path in the Brillouin zone identical as 
in Fig. [J] For N = 32, only one of the two lowest spin-wave 
branches [see second paragraph below Eq. (|5.9[) ] is shown. 



z spin components only on columns with j even. After 
this, not only compass- but also Heisenberg couplings 
on x bonds contribute to a^a^type terms in the LSW 



Hamiltonian. The resulting dispersion is: 

0JCz{k) = 4S\J (2\J Z \ + J kx + I kz f - (j kx + I kx f, 

(5.5) 

where 

{Jkx, ha] = {J x cosk x ,2Icosk a }. (5.6) 

Similarly, we derive the dispersion found by LSW ap- 
proximation for the lowest spin waves in Neel-like phases 
of the model, namely the G y and G z phase: 

uj Gy {k) = 4S'y / (ig + J kz - J kx f - (F % + J kx + J kz ) 2 , 

(5.7) 

u: Gz (k) = 4Sy/(2J z + U + J kx ) 2 - (J kx + I^) 2 . (5.8) 
Here we use as parameter definitions 

{Jkx, Jkz, I k } = {J x cosk x ,J z cosk z ,2I(cosk x + cosk z )}. 

(5.9) 

To illustrate these dispersions in the C z and G y phases, 
we show them along the path T— s-X— s-TU— ^F— 
in Fig. \TMb) (for I/J c = 0.2) and Fig. \M.c) (for 
I/J c = 0.6), respectively, with the same anisotropy pa- 
rameter <f> = 23/20 in both cases. In Fig. [TlTb) the cor- 
respondence between numerical results and the disper- 
sion Eq. (|5.5p is remarkable. Only in the energy range 
AE k > 5J C one notices a tiny discrepancy between nu- 
merical and LSW results. This can as in the F z phase be 
due to the interaction of singlc-magnon excitations with 
column flip excitations at energy cx L Z I. In Fig. H~4T c) 
the agreement between the LSW expression for the G y 
phase and the numerical results is also satisfactory, al- 
though finite-size effects are larger than in the case of 
Fig. [MTb) . The deviations from LSW theory for the G y 
phase are not due to excitations of energy oc L Z I (these 
are not properly defined in the G y phase) but are rather 
induced by the proximity of the C z — G y transition. 

In the phases corresponding to Figs. [14Tb) and WMc) 
and more generally in AF spin-gapped phases of the CH 
model, the lowest LSW branch is actually doubled due to 
the twofold ground state degeneracy in the classical limit: 
each branch contains the lowest spin waves above a lin- 
ear combination of classical ground states, which belongs 
to a given symmetry representation. In the example of 
the C z phase, this representation can be of momentum 
either k = T or k = X, and in the latter case the asso- 
ciated LSW dispersion is u>c z {k + X) instead of u)cz(k). 
A similar branch doubling occurs in Neel-like phases, re- 
sulting for the G a phase in a second branch of dispersion 
WGa(k + M) along with that of dispersion LOGa{k). 

In FM phases, branches are also doubled as a conse- 
quence of the Z2 symmetry of the ground state. This dou- 
bling does not appear on dispersion plots like Fig. [14Ta) 
since the two LSW branches, above ground states of iden- 
tical momentum Y, have the same dispersion; neverthe- 
less, and similarly as in AF phases, each spin-wave exci- 
tation found in exact diagonalization can be attributed 



16 



G 





N=8 




- H N=16 


0- 


N=18 


+- 


N=24 




-» N=32 


+ ■ 


+ LSW 





;.»=>- \ 


o-o(|)=n/4 
«3 ■<] (|)=7r/5 




a. 



0.05 1/fvJ 0.1 



0.2 



0.4 j /j 0.6 



0.8 



FIG. 15: (Color online) Anisotropy gap E a /J c of spin waves 
in the G z -phase as function of increasing J x / J z — tan(^), for 
Heisenberg amplitude / = J c /5 and various cluster sizes N. 
Inset: Size-scaling for the gap for: (f> — n/4 (isotropic case) 
and for <j> = tt/5 (finite gap). 



to one of these branches, thanks to its parity even or odd 
under time-reversal symmetry (or equivalcntly, in a phase 
with a easy axis, under the symmetry cr2 — > —o?). We 
plot both LSW branches in these figures; note also that 
one can deduce which translational symmetry is broken 
in the TL from the relative position in momentum space 
of both branches [e.g. the shift of q = X between both 
branches evidences a breaking of translation symmetry 
by e x in Fig.EJb)]. 

In an ordered phase of the model, the minimum of the 
dispersion of the lowest spin wave branch — or branches 
if taking into account the branch doubling in AF phases 
— is important for two reasons. First, the correspond- 
ing excitation energy (or spin gap) is to be compared to 
the energy oc L Z I of excitations mentioned above, which 
require a more detailed description given in the next Sec- 
tion. Second, when varying parameters of the model this 
spin gap vanishes at transitions with other phases, pro- 
vided the transition does not occur on the compass line 
/ = 0. A good example is the case of the G x — G z tran- 
sition, when J x is varied while / and J z are kept fixed 
and finite. In the G z phase the gap to spin waves E a 
is, for finite clusters, the lowest excitation energy in the 
sector of odd 55, found in representations of either 
k = r or k = M. It is shown in Fig. [T5] along with the 
corresponding LSW prediction from Eq. (|5.8[) . in func- 
tion of J x /J z . Even though finite-size effects arc not 
negligible away from the transition, i.e., for J x < J z , at- 
tempted scalings clearly indicate finite spin gap values in 
the TL which is comparable to the LSW theory result. 
Instead, at J x = J z , such a scaling confirms that there 
the spin gap vanishes in the TL. Within the LSW the- 
ory one finds that the spin gap vanishes at the transition 
as Cy/\J Z — J x \. The symmetry relation connecting each 
point of the G z phase, in the phase diagram Fig. [51 to a 



point of the G x phase and vice-versa, implies a relation 
between spin-wave dispersions in both phases; the LSW 
result for the G x phase is given by inverting, in Eq. (|5.8j) . 
first J x and J z , and second k x and k z (these dispersions 
being even functions of k x and k z ). 

One finds similar mode softening at the transition be- 
tween C z and G y phases, although these phases are not 
related to each other by any exact mapping. Here, if 
one approaches the transition from the side of the C z 
phase, by increasing \I/J Z \ with fixed J x /J z € (0 : 1), 
the dispersion wcz(fc) has a minimum at k = Y ac- 
cording to Eq. (|5.5|) [see Fig. iLlT b)] which vanishes for 
I/J z = —1/2. Similarly, in the G y phase and close 
enough to the transition towards the C z phase, the dis- 
persion given by Eq. (|5.7[) is characterized by a minimum 
at k = X (the other branch u>cz{k + M) has a minimum 
of equal value at k = Y) which vanishes at / = — J z /2 as 
well. Note that the closeness to one another of spin gap 
values E a = mine u)(k) in the three cases shown in Fig. PHI 
is accidental and is merely a consequence of the choice of 
the I/J z value for each case. The cases of C' z — F y and 
C' z — F x transitions, addressed in the previous Section, 
are characterized by similar mode softening, qualitatively 
reproduced within the LSW theory^ in each phase the 
spin gap corresponds (up to finite-size effects) to the low- 
est excitation energy seen in Figs. QT] and [12] at k = Y 
(F y or F x phases), or at k = T (C z phase). Interestingly, 
in an ordered phase $ but sufficiently close in the phase 
diagram to a transition line (except the line at / = 0) 
to another phase, the momentum ko corresponding to 
the minimum in the LSW dispersion w$(fc) allows one to 
deduce which translation symmetries are spontaneously 
broken in this other phase. 



VI. COLUMN-FLIP EXCITATIONS ON 
NANOCLUSTERS 

In the CH-model there is another distinct set of ele- 
mentary excitations, i.e., in addition to the spin- waves. 
These are the column-flip excitations that correspond to 
a reversal of all spins in a column of strongly coupled 
spins in the case \J Z \ > \J X \ and small Heisenberg am- 
plitude These excitations emerge from the macro- 
scopic ground state degeneracy of the original compass 
model and reflect the twofold ground state degeneracy of 
ordered phases selected by the Heisenberg interactions. 
Whereas spin-wave excitations yield essentially the same 
energy of order O(l) for small clusters (within exact diag- 
onalization) and in the TL, this is distinct for the column- 
flip excitations whose energy scales with a linear dimen- 
sion of the system. In the following paragraph we will an- 
alyze the column-flip excitations by means of an effective 
pseudospin model which we will derive from the original 
CH model using high-order perturbation theory. Then 
we will employ this effective model and will show how a 
quantum computation scheme involving the column-flip 
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excitations can be conceived. This requires the fulfilment 
of certain conditions on the low-energy excitation spec- 
trum, which we will eventually examine. 

A. Derivation of an effective columnar model 

We consider here finite clusters of size L x x L z with 
anisotropic CH interactions, assuming \J Z \ > \J X \ with- 
out loss of generality The amplitude / of Heiscnbcrg 
interactions is chosen finite but small compared to \J Z \, 
which is known (see Sec. 1111)) to be sufficient to lift 
the quasi-degeneracy between 2 ix columnar states. This 
splitting implies that, among those states, the (2 ix — 2) 
ones which do not correspond to the ground states of the 
selected ordered phase acquire finite excitation energies 
due to the finite value of /. Figure ITST a) shows energies 
of lowest eigenstates as function of I for a square clus- 
ter with edge length L x = L y = 4. For / small enough 
(< 0.15) a group of 16 distinct eigenstates lies below 
the lowest spin wave excitation — these 16 states have 
energies varying roughly linearly with /, with different 
branches corresponding to different slopes dE/dl. This 
feature is not specific to the FM CH model ( J z < 0); for 
the AF case (J z > 0) one sees in Fig. [TBTb) the same 
type of excitation branches with linear dependence of ex- 
citation energies on I when this quantity is small enough. 
The main difference is the somewhat larger splitting of 
excitation energies within a multiplet-branch in the AF 
case. 

We have seen in Sec. [IT] that an effective Hamiltonian 
-^coi ai l° ws valuable insight into the QCM (/ = 0) for 
the anisotropic case \J Z \ > \J X \. This effective model 
uses a formalism of pseudospins fj describing the colum- 
nar states forming the low-energy subspace. Here, to de- 
scribe the peculiar properties of the low-energy spectrum 
in the case where both / and J x are finite but small w.r.t. 
| J z | , we will derive a more general effective Hamiltonian 
H co i, expressed in the same pscudospin formalism. In 
the derivation we consider compass couplings on x bonds 
and Heisenberg couplings as perturbations. In the FM 
case (J z < 0), the resulting effective model is a ID XYZ 
Hamiltonian in terms of pseudospins Tj, 

H™ = -N\J Z \+ Yl E CTf-l+i' 

j=l...L x a£x,y,z 

where the coupling constants are given by: 

C z = L Z I, (6.2) 

Here J co i is given by Eq. (|2.9[) and depends again on clus- 
ter size and boundary conditions. The constant — N\ J z \ 
is the ground state energy of the unperturbed Hamilto- 
nian. This effective columnar Hamiltonian describes the 
structural and (quantum) dynamic properties of the low- 
energy, column-ordered states in CH nanoclusters. 



One can qualitatively interpret the difference between 
this Hamiltonian and -ff^oi obtained previously for 1 = 0, 
by listing the various roles played by Heisenberg cou- 
plings in the perturbation theory: 

(a) Most important are the Ia z a z couplings on hori- 
zontal bonds — they split the degeneracy of colum- 
nar states at first order in perturbation theory and 
contribute to the terms ex L z It z t z +1 which account 
for the ordering in the TL discussed in Sec. IIII1 

(b) The couplings Ia z a z on vertical bonds, instead, do 
not distinguish between the columnar states, but 
they contribute to their energy, either by a quantity 
L X LJ (with PBC) or L X (L Z - 1)/ (with OBC). 

(c) The transverse components 2I(a + a~ + c.c.) of 
Heisenberg couplings on horizontal bonds have to 
be added to terms J x o x o x when evaluating the 
transverse coupling amplitudes C x and C y . Here, 
not only t x t x terms but also (smaller) t v t v terms 
appear at order L z in perturbation theory. 

(d) Eventually, the transverse Heisenberg couplings 
H{p + o~ + c.c.) on vertical bonds have to be con- 
sidered a priori in the perturbative approach. For 
J z < they can be left out, since columnar 
states have spins ferromagnetically aligned within 
columns; but for J z > these couplings contribute 
to the columnar states As a result the effective 
Hamiltonian HQ^ is formally the sum of H™ and 
of an additive term: 

ffcoi = -'cd]>>;, (6.4) 

j 

I col = L z 2W2 7 (') /2 j z *' t (6 . 5) 

accounting for a single column flip. It appears at 
order L z /2 in perturbation. 

The strength \C y \ of the t»t* coupling is, for < |/| < 
\J X \, \J Z \, much smaller than that \C X \ of the t x t x cou- 
pling; the former vanishes for / — > 0, where one recovers 
the effective Hamiltonian i^oi- Both coupling strengths 
are, for L z large enough, much smaller than the one \C Z \ 
of the t z t z coupling. This allows for the latter coupling 
to split columnar states, so that the low-energy spectrum 
has the branch- like structure seen in Fig. \W[c) . For this 
system size there are three branches corresponding to 
the twofold degenerate ground state, the central one and 
the upper multiplct branch. These energy splittings are 
given by the number of domain walls, i.e., 0, 2, or 4 in 
the pseudospin chain. 

We comment here briefly to the implications of this ef- 
fective model for the interpretation of finite size data in 
previous Sections, in situations where \I\ <?C \ J X \ < \ J Z \. 
There, we stated that the reason why ordered phases are 
favored even with infinitesimal Heisenberg couplings ap- 
pears clearly with this effective description. And indeed, 
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FIG. 16: (Color online) Low-energy spectra of columnar ex- 
citations E/J c obtained by exact diagonalization of a L x = 
L y = 4 periodic cluster, as function of I / J c , for: (a) FM inter- 
actions (0 = 237r/20), and (b) AF interactions (<f) — 37r/20), 
respectively. Corresponding results obtained from the effec- 
tive Hamiltonian H™^ AF + H™^ AF for the same cluster and 
interaction parameters are shown for the FM (c) and AF case 
(d), respectively. Momenta of various eigenstates are indi- 
cated by symbols; continuous and dashed lines indicate even 
and odd states w.r.t. time reversal. 



from Eqs. (|6.3[) and (|6.5[) it is clear that while parameters 
Qx/y anc } / col vanish exponentially in the TL, C z does not 
(and even diverges in this limit). More quantitatively, for 
finite clusters one can expect that the order favored by 
effective t z t z couplings occurs when \C Z \ > \C X \ and 
\C Z \ > | /col |- Here the example shown in Fig. SI for 
I / J c = 10~ 6 , <fi = 7r/10, and L x = L z = 6, is instructive: 
the corresponding effective coupling amplitudes are given 
by I co i/(L z J z ) ~ 1.0 x lfr 18 , J col /(L z J z ) ~ 1.4 x 10" 6 , 
and C Z /(L Z J Z ) ~ 1.0 x 10~ 6 . Thus \C X \ is somewhat 
larger than C z , explaining the broad distribution of S z (k) 
over all momenta such that k z — tt, as on the compass 



line; but even for such small Hciscnberg amplitude /, C z 
is found not to be negligible relative to \C X \ and in con- 
sequence S z (k) is much larger at k = M than at other 
wave vectors. Similarly, the sharp peaks at I = in 
the fidelity plots of Fig. [TU] indicate that the Heisenberg 
coupling strength necessary for the C' z or G z order (de- 
pending on the sign of /) to develop on clusters consid- 
ered is even smaller than the resolution 0.005 J c chosen in 
that plot. According to the criterion C z 3> C x (with I co \ 
negligible in this regime) for <j> = — 37r/20 even a value 
| J | > 10 -3 J c is sufficient to develop long-range order. 

Until now we only considered the effective couplings 
found in leading order in perturbation for each pseu- 
dospin component x, y, z, and explaining the overall 
structure of the low-energy spectra in Figs. [ToTa) and 
I16f b). For a more accurate description of those, we in- 
clude in the effective Hamiltonian, besides H™^ AF , sub- 
leading terms H F ^^ AF found at second order in pertur- 
bation theory. The latter account for processes where two 
spins are flipped on nearest neighbor bonds. For J z < 
only horizontal bonds have to be considered; since the 
amplitude of these terms, for 1^0, depends on whether 
the pscudospins involved are parallel (amplitude J x ) or 
antiparallcl (amplitude J x + 21), the resulting nearest 
neighbor coupling is of the t z t z type. The part of the 
effective Hamiltonian stemming from these processes is 
(here for J z < 0): 
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f+i)- (6-6) 



The spectra in Figs. ITSf c) andfTfjld) correspond to effec- 
tive Hamiltonians #™+#™' and + H^J' , respec- 
tively — in the latter case differs from H FM ^ by 
the presence of an additional constant — iV— , account- 
ing for fluctuations on vertical bonds. The inclusion of 

^cot F auows f° r a much better agreement with the 
original spectra of Figs. fTtJl a) and lTBT b). regarding abso- 
lute energies and their dependence on I than if diagonal- 
izing H F ^f~^ AF alone. For instance, one recovers the fact 
that the dependence on I of energies of central and up- 
per branches is not simply linear but contains a (small) 
quadratic contribution. 

This effective model gives not only a good estimate 
of lowest excitation energies, but also the correct quan- 
tum numbers for the corresponding eigenstates within 
each branch. The two states of the lowest branch in 
Figs. flUla) and [TuT c). which in the TL become the 
twofold degenerate ground states of the F z phase, obvi- 
ously both have momentum T. Similarly, the two eigen- 
states forming the highest branch are linear combina- 
tions of states with a C z -like pattern, one at q = T 
and the other at q = X. The remaining, intermediate 
branch in Fig. I16f a) corresponds to the subspace gener- 
ated by 2 4 — 4 = 12 columnar states such that, in terms 
of pseudospins, J2i T i T i+i = ; tms branch contains 
I T1T2I3I4), I tit2T3|4), and I tiMsW, plus for each 
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of those the three states obtained by translations along 
the pseudospin chain. The effective model allows to un- 
derstand why three eigenstates of this branch are found 
in each representation of momentum such that q z = 0. 
Also the dispersion within a branch is well reproduced by 
the effective model, however the splittings are too small 
for this dispersion to be visible in the spectra of Fig. [T5] 
In the AF CH case (for J z > 0), the presence in the 
effective Hamiltonian of H x ol has a significant impact 
on properties of the intermediate branch, clearly visi- 
ble in Figs. [TBlb) and IT6l' d): since the amplitude / co i 
is much larger than C x ' v for / large enough (/ > 0.05 J c 
in the present example), the energy dispersion within this 
branch is dominated by J co i in this interaction range, and 
much broader than in the FM case [Fig.[T5Ja,c)] for equal 
value of |/|. 

We make here two important remarks on the general- 
ization of features above to other clusters: 

(a) In the previously discussed case of a periodic clus- 
ter with L x = 4, only one intermediate branch with 
a large quasi-degeneracy is found. This is spe- 
cific to this case, and for clusters with larger L x 
one has several such branches — for instance with 
L x = 6 and PBC one has two of them; each one is 
generated by 30 columnar states where J2i T i T i+i 
takes the value —2 or +2 respectively. The maxi- 
mal number of states per branch, which is given by 
2C£ with p (one of) the even integer(s) closest to 
L x /2, grows thus exponentially with the number of 
columns. 

(b) The obtained branch structure of the low-energy 
spectrum is a property of finite and small enough 
clusters, but does not require periodic boundaries. 
We have verified that for open clusters the descrip- 
tion of these systems with perturbative techniques, 
as in the periodic case, is possible and efficient as 
well, but the obtained effective Hamiltonian has a 
more complicated form: amplitudes of transverse 
couplings t"t" +1 depend on whether one of the 
columns (j or j + 1) is or not at the border of 
the cluster. The amplitudes of dominant t z t z cou- 
plings are the same as for PBCs, but the number of 
branches split off by those is larger (L x — 1) since 
odd numbers of domain walls are allowed in con- 
trast to the periodic case. 

B. Quantum computing scheme based on 
quasi-degenerate columnar states 

We have seen that the low-energy spectrum of open 
L x x L z clusters formed by er-spins that interact via 
CH couplings can be well reproduced, for small enough 
Heisenberg amplitudes and sufficiently anisotropic com- 
pass couplings, by an effective XYZ-type model of pseu- 
dospins r = 1/2. A remarkable feature of this spec- 
trum is the subdivision of the 2 Lm columnar states, form- 



ing a low-energy subspace selected by dominant compass 
interactions, in multiplet branches of quasi-degenerate 
states. Some branches have a semi-macroscopic degen- 
eracy, while the lowest branch contains only the two de- 
generate ground states selected by Heisenberg perturba- 
tions. The ground states are characteristic for the or- 
der of the respective phase and the Z2-symmetry of the 
model. The effective model allows, thanks to the high 
anisotropy of its coefficients (\C Z \ ^> \C X \,\C V \), for a 
good understanding of the number and nature of excited 
states within a branch, which can be described in terms 
of pseudospins 1/2. 

This suggests that one can control the quantum state 
of the system and possibly realize elementary operations 
of quantum computing by using a subspace of quasi- 
degenerate states. The starting point is to excite the 
system into a given branch, by flipping given columns; 
then one can initialize the quantum computer by plac- 
ing the system in a state where some pseudospins (the 
qubits of the computer) are highly entangled; and, af- 
ter this, perform quantum operations by acting on these 
qubits. For concreteness, we specify a system of inter- 
est: a rectangular, open cluster of L x x L z spins (we set 
in our example for concreteness L x = 5 for the rest of 
this paragraph). Concerning interaction parameters, we 
choose FM compass (J z < J z < 0) and AF Heisenberg 
(/ > 0) couplings. This choice is, within the framework 
of the CH model, the closest one to possible realizations 
using arrays of (pseudo) spins coupled by dipolar interac- 
tions [see Eq. (|1.1[) ]. One of the two lowest eigenstates 
of the system (these are quasi-degenerate ground states) 
can be expressed, as in Fig. ll7f a), in terms of pseudospins 
n by: 

ICS) = | ti+ataitts). (6.7) 

By flipping the j = 3 column, i.e., the pseudospin 73, one 
obtains the state shown in Fig. [TWb) : 

|EX) = I tiUhUW- (6.8) 

This state belongs, in a limit where C x ' v < C z , to 
a branch of 12 quasi-degenerate states; a basis of this 
branch consists of eigenstates of t z operators, such that 
t* = if = ±1 and for exactly two values j 6 [1; 4] one has 
ijij+i = — 1. We will consider from now on operations 
within this subspace, which we call central branch since 
the number of domain walls is half of the maximum al- 
lowed (L x — 1). Actually wc will further restrict the work 
subspace, i.e. the subspace spanned by all possible states 
available via operations on qubits, by imposing that the 
pseudospins f\ , 13 and T5 remain in the same state as in 
I EX). This leaves 2 pseudospins t 2 and 14 which can be 
operated by pulses acting coherently on all spins of the 
corresponding column, while keeping the system in this 
branch of states. One can first initialize this 2-qubit com- 
puter by applying Hadamard gates on pseudospins t 2 and 
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FIG. 17: Pseudospins of a CH nanocluster representing L x — 
5 columnar spin chains in the case | J z \ S> | J x \, Panels (a)- 
(d) show the successive steps in the initialization and utiliza- 
tion of two essentially decoupled protected qubits: (a) the sys- 
tem is in one of the two quasi-degenerate ground states, |GS), 
such that for each pseudospin j one has rj = (— 1) 3 ~ ; (b) af- 
ter flipping the column (or pseudospin) j — 3, the system is 
in an excited state |EX) belonging to the central branch of 
columnar excitations; (c) the qubits encoded by pseudospins 
T2 and are initialized in the state HzHilEX) , where Hj ro- 
tates the pseudospin fj along the y axis in pseudospin space; 
(d) example of a state produced from the former by an oper- 
ation on pseudospin fb. 



T4 , resulting in the following state [shown on Fig. ITTTc)] : 

H 2 H 4 |EX >= |0 >= ^ E ltiT 2 z | 3 T 4 z t5>, (6.9) 
T|iT|et;-l 

with Hj corresponding to the operator (tJ+tJ ) / v2- Fur- 
ther actions on either of the columns j ' = 2 and j = 4 
(such that the system remains in the work subspace) will 
correspond to elementary operations of the corresponding 
qubits, leading to intermediate states as the one repre- 
sented on Fig. flTTd). This scheme can be extended to 
systems with more qubits, namely p when using an open 
cluster with L x = 2p + 1 columns; in this more general 
case I EX) is an eigenstate of all rj operators such that 
T j T j+2 = — 1 f° r au 0( id columns j = 2k — 1 (k G 
qubits arc then encoded in the remaining columns of in- 
dex j even. 

This encoding scheme presents several advantages: 
first the qubits are semi- locally defined, i.e., each qubit 
corresponds to degrees of freedom of a specific column 
and can be manipulated by a field acting on this column. 
Second, the structure of columnar excited states makes 
them robust against local noise. When considering local 
perturbations to the CH Hamiltonian, for instance of the 
form randomly distributed (and small) j22, these can re- 
verse columns only at order L z in perturbation theory. 
Qubits based on columnar excitations are thus intrin- 
sically, as in the case of topological quantum comput- 
ing, much more robust to such perturbations than qubits 



which would be defined locally, e.g. on a single spin. Be- 
sides, the choice of using the work subspace described 
above has a further advantage. To see this, we consider 
the operator Qj that flips the whole column j = 2k where 
the fcth qubit is encoded; the commutator of this opera- 
tor with rl is given by Eq. (|2.7[) , from which we see that 
within the work subspace the a = z term of the sum van- 
ishes (this is because in this subspace ofj_i + = 
on all rows i). In consequence, expectation values of 
[Qj , %] are expected to be smaller within this subspace 
than within any other quasi-degenerate branch of colum- 
nar states. Thus the scheme described here appears as 
the optimal way to encode qubits using low-energy eigen- 
states of CH nanoclusters. 



C. Protection against relaxation and decoherence 

The quantum computation scheme presented above 
would work perfectly if: (i) within the branch involved in 
the scheme, eigenstates were exactly degenerate, and (ii) 
the quantum gates realizing the transition from one state 
to another could be implemented perfectly by physical 
operations. The latter condition means that one could 
apply a spatially-resolved magnetic field, focused on a 
given column of spins, for a fixed time t\; and, noting 
the initial state \A) (belonging to the work subspace) at 
t = 0, the state e l to 1 #(*)<** | A) obtained after this opera- 
tion would still belong to this subspace. Yet, in a physical 
system with CH-typc interactions, such operations would 
necessarily involve relaxation (energy can be exchanged 
with the environment and via the magnetic field so that 
the system relaxes to its ground state) and decoherence. 
The relaxation rate(s) in the system depend(s) on its 
details, but an important factor on which one can focus 
within the CH model is the energy difference between the 
branch of interest and the more conventional single-spin 
excitations, i.e., spin waves. Indeed, when one applies 
an imperfect field with the aim to flip a whole column, 
one spin of the column can remain unflipped, or a spin 
of a neighboring column can be unintendedly flipped; in 
both cases this schematically corresponds to single-spin 
excitations, as the spin waves seen in Sec. [Vj But if the 
gap to spin waves is large enough, these processes should 
have a negligible impact on the evolution of the system at 
moderate time scales. In consequence, we will impose, as 
a criterion for the robustness of such a scheme to such re- 
laxation processes, that the lowest single-spin excitation 
has an energy E a higher (if possible, much higher) than 
the energy of the highest eigenstate in the central 
branch, containing the work subspace. 

We examine the dependence of both excitation ener- 
gies first on interaction parameters, and second on cluster 
dimensions L x and L z . Figure [TSt a) shows the depen- 
dence of (and of E~ , excitation energy to the lowest 
eigenstate in the same branch) in function of I j J c for 
cluster parameters [L x — 3, L z — 4) and two distinct 
values of 4>. In these examples, even for the largest value 
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FIG. 18: (Color online) Lowest single spin-flip excitation en- 
ergy E a in the C z -phase in comparison with the maximal 
(E^~) and minimal (E^) energies of the central column- flip 
excitation branch [see text] for open clusters: (a) as function 
of I/Jc, for fixed L x = 3, L z = 4, and either <f> = ll7r/10 or 
<f> = 237r/20 — here and E^ are undistinguishable; (b) as 
function of L z , for fixed I / J c = 0.1, cj> — llfr/lO, and either 
L x = 3 or L x — 5. 



of J x /J z (i.e., 4> ~ 237r/20) the dispersion within this 
branch is negligible with respect to its lowest excitation 
energy (\E+ — E~\ <C E~); the linear dependence of Ef 
on / and the weak dependence on <p are consistent with 
the fact that these are columnar excitations, described 
by the pseudospin formalism of Sec. IVI Al In contrast 
the gap E a to single-spin excitations is much more sen- 
sitive to J x /J 2 ; it is reduced when interactions become 
less anisotropic, as in the bulk case discussed in Sec. |V] 
- and has only weak dependancc on I / J c in the regime 
considered. We stress that E a is, as in Sec. |V] the lowest 
single-spin excitation, but here due to open boundaries 
it is roughly twice as small as the value expected, with 
identical interaction parameters, from Eq. (|5 . 5|) . Indeed, 
open boundaries allow for edge modes of lower energy 
than the bulk spin wave modes, since flipping a spin at 
one (horizontal) edge frustrates only one z-bond compass 
coupling instead of two. 



The dependence of column-flip excitations on cluster 
dimensions is displayed in Fig. fT5Tb). The figure shows 
for fixed interaction parameters the values of E+, E~, 
and E a as function of L z for two values L x = 3 and 5 
compatible with the scheme described in Sec. IVI Bl The 
linearity of E^r in L z , expected from perturbation the- 
ory when assuming l/),),/^ <C |J Z |, is clearly verified. 

The slope is about twice as large for L x = 5 as 

for L x = 3. Indeed, in the first case the central branch 
is built on columnar states which contain two domain 
walls (i.e., in which t*t* +1 = 1 for exactly two values 
of i = 1,..,L X — 1), while in the latter case the corre- 
sponding branch is built on columnar states with only one 
domain wall. In contrast, we see that E a is roughly size- 
independent (only for L z = 2 its value differs noticeably 
from those for longer columns), and reduced of ~ 25% 
from the value 2\J Z \ expected in the \I\, | J x \ <C \J Z \ limit 
- an adaptation of LSW theory to finite clusters could 
provide a closer estimation taking into account the effect 
of Heisenberg and x-compass couplings. 

Based on these features, we define a column- flip regime 
where the columnar excitations corresponding to the cen- 
tral branch have lower energy than the lowest single-spin 
excitation, i.e., < E a ; the complementary case is 
called spin-wave regime. In Fig. I18f a) we see that the 
extent of the column-flip regime, in terms of I / J c , is de- 
creased with increasing J x / J z for (L x , L z ) = (3, 4) — this 
regime corresponds to I / J c < 0.20(1) for <f> = Htt/10 and 
to I / J c < 0.17(1) for <j> = 23-7I-/20. Increasing cluster di- 
mensions reduces the extent of the column-flip regime in 
parameter space, due to the strong dependence, seen in 
Fig. 1Mb), of E+ on both L x and L z parameters. 

Importantly, the column length has also a strong im- 
pact on E+ — E~: if e.g. L x = 3, the four states of this 
branch are quasi-degenerate (the splitting is not visible 
to the eye) for L z > 3, while they are split by energies 
of < 0.05 J c for L z = 2. The same trend is observed for 
L x = 5, but the higher number of states in that branch 
results in a significantly broader dispersion, ~ 0.2 J c for 
L z = 2 and decreasing rapidly with increasing column 
length but still visible in the figure for L z = 3, 4. This 
is in agreement with the perturbation theory approach 
in Sec. IVI Al according to which the splitting within a 
branch, governed by parameters C\ 
exponentially with increasing L z . 

We can now make use of these results and formulate 
the conditions for such a system to be in the column-flip 
regime, and simultaneously require that the dispersion 
E+ — E~ <C E~ , so that related decohcrcncc effects are 
as small as possible. The first condition requires rather 
short columns (because of the scaling E^r cx L Z I), while 
the second condition requires sufficiently long columns 
[L z > 2 with the interaction parameters chosen for 
Fig. [TST b)]. In the L x = 3 case, which might allow us 
to engineer a one-qubit device, for values L z = 3,4,5 
both conditions are fulfilled; but in the L x = 5 case, 
which would correspond to a two-qubit device, these con- 
ditions are fulfilled only for L z = 3. Of course, the scal- 
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ing oc L Z I implies that not only the column length 
but also the amplitude of Heisenberg couplings must be 
rather small; and from Sec. |V] we know that small values 
of J x /J z (he., strongly anisotropic compass couplings) 
are preferable, to keep these excitation energies below 
those of single-spin flips. Thus, in principle, one can en- 
gineer clusters with large dimensions, i.e., large number 
of qubits, with realizations of the CH model where J x / J z 
and I / J z are tunable at wish. 

One can consider these conditions with a specific phys- 
ical system in mind, namely arrays of quantum spins 
with dipolar interactions (for instance, representing NV 
centers as discussed in Sec. HJ; there, the choice of 
interaction- and size parameters is restricted further. If 
we neglect the role of interactions beyond nearest neigh- 
bor effective spins,— the Heisenberg- type couplings on 
vertical bonds and on horizontal bonds have distinct 
(negative) amplitudes, J z /3 and J x /3 respectively, in- 
stead of a uniform amplitude /; and the anisotropy ra- 
tio J x /J z is fixed by the geometry of the array to ( 3 , 
with £ = c/a 6 (0; 1) the ratio between vertical (c) and 
horizontal (a) bond lengths. Note that £ is the only pa- 
rameter governing all ratios between the various coupling 
amplitudes. 

One can imagine an array of spins with dipolar cou- 
plings, restricted by some screening mechanism to nearest 
neighbors; and such that, as in Fig.fTSTb), I x /J x ~ —0.33 
(with I x the Heiscnbcrg-type amplitude on horizontal 
bonds) and J x /J z ~ 0.33. The corresponding aspect 
ratio would be £ = 0.69, but the main difference between 
this situation and the case of Fig. flST b) would reside 
in the amplitude of Heisenberg-type couplings on verti- 
cal bonds (I z ~ 0.3J C in the dipolar case compared to 
I = I z = I x ~ 0.1J C in the CH case). This difference 
should not play a significant role on values of E^ since 
the corresponding states of this branch have ferromag- 
netically aligned columns; but in the dipolar case with 
C = 0.69, we checked that the gap to lowest single-spin 
excitations would then be reduced compared to E a seen 
in Fig. H5T b). by a factor of ~ 1.5 only, and branches 
of quasi-degenerate excitations would still be present in 
the low-energy spectrum of e.g. the (L x = 3, L z = 4) 
cluster. Thus, in principle, an array with this or smaller 
values of £ would allow to define a one- or two-qubit sys- 
tem, assuming that such highly anisotropic arrays, with 
dominant dipolar interactions, are conceivable on an ex- 
perimental point of view. 



VII. SUMMARY AND DISCUSSION 

In this work we have investigated the ground states 
and elementary excitations of the compass-Heisenberg 
model formed by quantum spins. The compass model 
is characterized by a macroscopic ground state degener- 
acy, therefore one of our central goals is to explore what 
happens with this huge degeneracy when the system is 
not perfect but is exposed to other perturbativc interac- 



tions, which we assume here to be of Heisenberg type. 
Compass interactions can arise in the description of a 
variety of strongly correlated electronic systems, ranging 
from orbitally degenerate Mott insulators to cold atoms 
or ions in optical lattices, the latter having attracted at- 
tention in the quest for possible realizations of quantum 
computing devices. Besides, we pointed out that the 
compass-Heisenberg model can be seen as a short-ranged 
version of a Hamiltonian for NV centers coupled by dipo- 
lar interactions, and also studied intensively recently with 
motivations from quantum information. 

We first analyzed the zero-temperature phase diagram 
of this model, using analytical approaches and numerical 
exact diagonalization of the Hamiltonian. We have found 
that a feature characteristic of the compass model, the 
semi-macroscopic ground state degeneracy in the ther- 
modynamic limit (TL) is lifted in presence of Heisenberg 
interactions even for infinitesimally small amplitude I of 
the Heisenberg coupling. As a result, the phase diagram 
contains various ordered phases, either ferromagnetic or 
antiferromagnetic. The latter include different columnar 
ordered phases. Due to the anisotropy of interactions in 
spin space, i.e., except for special values of the parame- 
ters 1 1 J z and J x j J z , these phases have an easy spin axis, 
and the ground state degeneracy in the TL is twofold. 
Transitions between these phases occur for coupling am- 
plitudes either equal or very close to the corresponding 
values in the classical (large spin) limit of the model — 
in the first case, this apparent insensitivity to quantum 
fluctuations is a consequence of extra symmetries of the 
Hamiltonian at particular transition lines. In the second 
case, we presented a rather precise perturbative evalua- 
tion of the quantum corrections to the phase boundaries, 
i.e., consistent with the shifts obtained by exact diago- 
nalization. 

The phase transitions in the compass-Heisenberg 
model are continuous on the finite systems studied; ten- 
tative size scalings indicate that they may become of first 
order in the TL, but alternatively they may keep a con- 
tinuous character, since they are characterized by the 
presence of gapless excitations in contrast to the ordered 
phases selected elsewhere. The modes becoming soft at 
the transitions can be of two types: for those occurring at 
finite / these are spin waves, for which the dispersion in 
ordered phases can be well described in the linear approx- 
imation. The case of transitions at / — between two 
magnetically ordered phases is specific: such a transition 
corresponds to a level crossing between a multitude of 
columnar states, which allow to define excitations char- 
acteristic of the compass-Heisenberg model. 

These column-flip excitations consist of flips of all 
spins within a column, assuming that the dominant 
compass couplings are those on vertical bonds. The 
corresponding excited states, which belong to the sub- 
space spanned by ground states of the compass model 
at I = 0, are for small but finite / split off proportion- 
ally to this amplitude and to the size of columns. This 
also implies they are pushed up to high energies in the 
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TL. In small nanoscale systems, however, column flips 
can be the lowest excitations, i.e., they may lie inside 
the anisotropy gap of spin waves. Column-flip excita- 
tions also have the remarkable property of being grouped 
into multiplet-branchcs of quasi-degenerate states — the 
splitting within a branch decreases exponentially with 
increasing column length and the number of states in 
certain branches grows exponentially with the number of 
columns. These features, and more generally the proper- 
ties of these excitations, are described in detail by an ef- 
fective one-dimensional model, which we derived in high- 
order perturbation theory. This effective model takes the 
form of a XYZ-type Hamiltonian between nearest neigh- 
bor pseudospins t = 1/2. 

The effective low-energy Hamiltonian allows, in the 
regime where this perturbation theory applies, to de- 
scribe the quantum dynamics of a finite cluster of L x x L z 
spins, with transitions from the ground state to excited 
columnar states and between the latter ones. Based on 
this, we propose a novel type of quantum computing 
device, where qubits are physically encoded in specific 
columns of a cluster, and where the work subspace is 
formed by a quasi-degenerate excitation branch. The 
encoding in columns instead of single spins renders the 
system fault tolerant, similar to proposals for topological 
quantum computing — although a non-trivial topology is 
not required here. We suggest here a possible realization 
of this type of encoding by means of rectangular arrays 
of quantum spins (each spin representing e.g. a NV cen- 
ter in a diamond matrix) coupled by dipolar interactions 
rapidly decaying with distance. The truncation of these 
interactions to first neighbors represents, up to minor de- 
tails, a particular case of the compass-Heisenberg model 
with ferromagnetic compass couplings. Thus the pecu- 
liar features of low-energy excitations in the latter model 
could also be encountered in this more realistic system. 

The ferromagnetic nature of dominant compass cou- 
plings presents an advantage for such a realization, as 
it allows one for an easy manipulation of columns by 
an external field. A device would have to fulfil several 
conditions on the multiplct-branch containing the work 
subspace: (i) it should be well separated from single- 
spin excitations, and (ii) the energy splitting within this 
branch should be sufficiently small. We analyzed these 
conditions, finding that they impose restrictions on the 
dimensions of the array. To satisfy both conditions, an 



optimum has to be found for the length of columns, while 
the number of columns, determining the number of hy- 
pothctically available qubits, has to be limited to satisfy 
the former condition. This number can nevertheless be 
increased by varying the geometry (namely the aspect 
ratio £) of the array, although an unrealistic aspect ratio 
corresponding to quasi-decoupled columns may lead to 
other (dccohercncc-relatcd) problems. 

A possible development of this study would be to sim- 
ulate the time evolution of arrays of spins within the 
compass-Heisenberg model and to estimate decoherence 
and relaxation times. This involves to compute the re- 
duced density matrix corresponding to the pseudospins 
defining the qubits of the system, while other degrees of 
freedom are traced out. The time evolution can be stud- 
ied in the framework of the compass-Heisenberg Hamilto- 
nian itself, where one can also add perturbing terms ac- 
counting for unavoidable noise effects, and it is also pos- 
sible to model elementary operations by including time- 
dependent fields centered on specific columns. In the lat- 
ter context the effective Hamiltonian represents a great 
advantage as it allows to simulate the time-evolution of 
qubits much more effectively, i.e., compared to the origi- 
nal compass-Heisenberg model. 

Another direction to follow, more closely connected 
to arrays of NV centers, would be to consider, instead 
of the (short-ranged) CH interactions, the (power-law- 
decaying) dipolar interactions and reexamine the condi- 
tions for a similar encoding. Eventually, taking a more 
theoretical point of view, the apparently contradicting 
features seen at the quantum phase transitions occur- 
ring at / 7^ (mode softening in the spin-wave spec- 
trum, versus size scalings indicating a conventional first 
order scenario in the TL) call for further work using com- 
plementary approaches, especially for the transitions not 
characterized by additional symmetry on the transition 
line. 
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